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^SJ , Abstract: We compute Yukawa couplings involving chiral matter fields in toroidal com- 

\ pactifications of higher dimensional super- Yang-Mills theory with magnetic fluxes. Specif- 

■r^ij- I ically we focus on toroidal compactifications of D=10 super- Yang-Mills theory, which may 

' be obtained as the low-energy limit of Type I, Type II or Heterotic strings. Chirality is ob- 

tained by turning on constant magnetic fluxes in each of the 2-tori. Our results are general 
and may as well be applied to lower D = 6, 8 dimensional field theories. We solve Dirac and 
D . Laplace equations to find out the explicit form of wavefunctions in extra dimensions. The 

Yukawa couplings are computed as overlap integrals of two Weyl fermions and one complex 
scalar over the compact dimensions. In the case of Type IIB (or Type I) string theories, the 
^ ' models are T-dual to (orientifolded) Type IIA with D6-branes intersecting at angles. These 

d ' theories may have phenomenological relevance since particular models with SM group and 

three quark-lepton generations have been recently constructed. We find that the Yukawa 
couplings so obtained are described by Riemann i?-functions, which depend on the complex 
structure and Wilson line backgrounds. Different patterns of Yukawa textures are possible 
depending on the values of these backgrounds. We discuss the matching of these results 
with the analogous computation in models with intersecting D6-branes. Whereas in the 
latter case a string computation is required, in our case only field theory is needed. 



Keywords: D-branes, Yukawa couplings. Extra Dimensions, String Phenomenology. 



Mirror Symmetry 



Contents 



i 



i 



Introduction 

Yukawa couplings in Kaluza-Klein theories 

!T| D = 10 N=l Super Yang- Mills compactifications with magnetic fluxes 
?2l Models with fluxes on D = 6, 8 dimensions 



Toroidal wavefunctions I: Abelian Wilson lines 

7\\ Eigenfunctions of the Dirac equation on 



3.1.l| Abelian gauge field 
3.1.21 Dirac zero modes 
Normalization 



3.1.3 



3.2 



Chiral matter eigenfunctions 
3.2.1 Fermions in bifundamentals 



3.3 



3.4 



Eigenfunctions of the Laplace equation 
Theta functions as wavefunctions 
Generalization to T^" 



Toroidal wavefunctions II: non- Abelian Wilson lines 
[4.1| Non-Abelian gauge groups 
4.2 Fermions in bifundamentals 



4.3 Eigenfunctions of the Laplace equation 

4.4 Generalization to T^" 
4.4.1 Factorizable tori 
4.4.21 General tori 



Computing Yukawa couplings 

H Computing Yukawas on a 
|5.1.1| Abelian Wilson lines 
|5.1.2| Non-Abelian Wilson lines 
Higher dimensional tori 
Yukawas in supersymmetric models 



5.2 



5.3 



D-branes of lower dimension 

Comparison with intersecting brane computations 
|7.1| Intersecting brane Yukawas 
Matching by T-duality 

Horizontal T-duality 



7.2 



7.2.1 



7.2.2 Tilted T-duality 



7.3 



7.4 



Chiral fields and T-duality 

The Kahler metrics of chiral fields revisited 



10 

H 

11 

12 
15 
16 
17 
19 

m 



26 

m 



37 

|8 
|9 
43 
44 
45 

46 

48 
|9 
1511 



- 1 - 



m. A 3-generation MSSM-like Model 



56 



|9|. Final comments and conclusions 

\A[ Dimensional Reduction of N=l Super Yang-Mills 

|B|. Fluxes and supersymmetry 



B.l Hermitian Yang-Mills Equations 



59 
61 
64 

H 



B.2 D-brane interpretation 



B.3 Mirror symmetry and string corrections 



1. Introduction 

One of the most outstanding puzzles of the standard model (SM) of particle physics is the 
structure of the Yukawa couplings between the Higgs field and the SM fermions. A correct 
description of the observed masses and mixing of quarks and leptons seems to require very 
different values for the Yukawa coupling constants for the different generations. Although 
many approaches have been attempted to describe the hierarchical structure of Yukawa 
couplings it is fair to say that we do not have at the moment a compelling theory for 
quark and lepton masses. Thus the search for theoretical schemes which could explain this 
Yukawa structure is of the utmost importance. 

In the last 25 years the idea that there could be more than four dimensions has been 
pursued intensively, particularly due the study of string theory which is naturally defined 
in 10 or 11 dimensions. Extra dimensions offer in principle the possibility of computing 
Yukawa couplings in terms of the extra-dimensional geography. Indeed, starting from a 
(D-|-4)-dimensional field theory and compactifying D dimensions one may get massless 
modes with factorized wavefunctions ip{x) x (p{y), with x{y) denoting Minkowski and ex- 
tra dimensions respectively. Gauge boson components A* in extra dimensions give rise to 
scalars at low energies and Yukawa couplings are thus expected to appear upon compactifi- 
cation from the higher dimensional gauge vertex interaction A^"^Tm^- Qualitatively the 
Yukawa coupling constants may be computed from overlap integrals over extra dimensions 
of the form [1,2] 

Yij = I dy^'MyHMMy) ■ (i-i) 

Unfortunately it is not easy to find compactifications in which: 

• One has chiral fermions, as required in the SM 

• One can solve the Dirac and Laplace equations to obtain explicit expressions for the 
wavefunctions in extra dimensions 

• One can explicitly work out the overlap integrals and actually compute the Yukawa 
couplings 



Our main goal in the present article is to perform such a type of computation in a 
class of theories of potential phenomenological interest. We consider as our starting point 
10-dimensional super- Yang-Mills (SYM) theory as the best motivated extra dimensional 
field theory, since it appears in the low-energy limit of Type I, Type IIB and heterotic 
string theories. However one can easily apply the results (i.e. the solutions of Dirac 
and Laplace equations as well as the computation of overlap integrals) to other lower 
dimensional D = 6,8, as we discuss below. 

We compactify D=10 SYM on a 6-torus and, in order to obtain chiral fermions, 
we add constant magnetic flux through the torus [3-6]. The fact that one may obtain 
chiral fermions in the presence of explicit gauge field backgrounds is well known [7] . How- 
ever, in the present case we consider a simple toroidal geometry with constant gauge field 
backgrounds, so that one can explicitly find the eigenfunctions of the Dirac and Laplace 
equations in extra dimensions. ^ We find that in the case of a factorized torus T^" the 
eigenfunctions are proportional to products of n Jacobi theta-functions with characteristics. 
The profile of the corresponding wavefunction densities in extra dimensions is Gaussian- 
like, with the location of the maximum controlled by the possible presence of Wilson-line 
backgrounds (see figures 3, 4 and 5). 

From these wavefunctions one can explicitly perform the overlap integrals and obtain 
the Yukawa couplings. The Yukawas so obtained are again products of n iJ-functions but 
this time they depend both on the complex structure moduli of the tori as well as on 
the possible Wilson line backgrounds present in the compactification. In the case of a 
non-factorizable 6-torus one obtains analogous results but now one has instead generalized 
Riemann theta functions. 

One of the main results of the present paper is the final expression for Yukawa cou- 
plings. In the case of a factorized 2n-torus, the Yukawa couplings between three fields 
labeled by integers i,j,k have the following general form: 



^ijk — 54+2 
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Here 54+271 is the (4 -|- 2n)-dimensional gauge coupling constant in the dilute fiux approx- 
imation^, T^^\ A^^^ are the complex structure and the areas of the r = 1,2,3 compact 
2-tori whereas C^^^ describe the dependence on the Wilson line variables. Yir ^"^^ 

(r) 

number of generations of particles labeled by i. The quantity 9j is the magnetic fiux felt 
by such particle in the r^^ 2-torus, which also signals the scale of new physics in this sector, 
more precisely the mass of the lightest massive replica of the chiral fermion i. H^^^ is a 
known function (see eq. ( 5.32D ) of Wilson lines, complex structure and magnetic fluxes. It 
vanishes in the absence of Wilson lines. Finally, "& are standard Jacobi theta functions with 



characteristics, with the argument 5^^^ = + j/Ij' + k/I^' , being the only 'flavour' 



(r) 



r(r) 



ijk 



^Toroidal compactifications with constant field strengiit are also known as 'Torons' [8-11]. The connec- 
tion with D-brane physics has been adressed in [12-16]. For other related computations see [17]. 

^This coincides with the gauge couplings of the different gauge groups up to flux dependent corrections 
which are suppressed in the large volume limit. 
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(i, j, k) dependence of the Yukawa coupling. In a SUSY compactification the (holomorphic) 
superpotential should be identified with this product of t?-functions, whereas the rest of 
the factors should come from D-term normalization of three-point functions. 

As we said, the natural setting for studying compactifications of D=10 super Yang- 
Mills is string theory, since it appears naturally in the effective field theory limit of Type 
I, Type JIB (with D9-branes) and heterotic strings. In fact, although the mentioned com- 
pactification and computations may be performed without any reference to string theory, 
in the present paper we will have always in mind the string theory point of view. In partic- 
ular it is well known that Type IIB theory with D9-branes compactified on a (magnetized) 
is T-dual to Type IIA theory with D6-branes wrapping 3-cycles on and intersecting 
at angles. ^ In recent years string models with intersecting D-branes have been actively 
pursued in order to obtain realistic models of particle physics [4,6, 18—20]. The main mo- 
tivation for this is that in these schemes chiral fermions naturally appear at the points in 
compact space in which Dp-branes intersect. Furthermore, since Dp-branes may intersect 
multiple times, this gives a rationale for the family replication present in the SM. A num- 
ber of semirealistic models with the massless fermion spectrum of the SM or some J\f = I 
extension with three generations have been obtained [20-25] ^. 

In fact one of the motivations for the present investigation was to check how T-duality 
maps the results obtained for intersecting D-branc models to the effective action of Type IIB 
with D9-branes and magnetic fluxes. In particular, in ref. [27] we computed the 'classical' 
contribution to Yukawa couplings in Type IIA models with intersecting D6-branes. In this 
case one has to perform a sum over string worldsheet instanton contributions to obtain the 
final expression of Yukawa couplings, a pure stringy (non field-theoretical) computation. 
On the other hand T-duality tells us that cqTiivalcnt Yukawa couplings should be obtained 
for the case of Type IIB theory compactified on with magnetic fluxes. In this second case 
the computation is purely field theoretical since one just starts from D = 10 SYM theory, 
compactifies and performs the overlap integral to obtain explicitly the Yukawa couplings, 
as described above. Thus on this flux side of the T-duality the computation is an exercise 
on Kaluza-Klein theory. Indeed we have found that the Yukawa couplings computed in 
both T-dual theories agree after appropriate transformations of the moduli and Wilson line 
variables. This is a non-trivial check given the completely different origin (stringy on the 
intersection side, field theoretical on the flux side) of both computations. Furthermore, 
the computation obtained in the present paper goes beyond the results obtained in [27] 
since we are also able to obtain all normalization factors and coefficients appearing in the 
Yukawa coupling. In the intersection side in order to obtain the normalization factors one 
has to compute the quantum piece of the appropriate string correlators as in refs. [28,29]. 
We find agreement with those normalization factors in the limit of small intersecting angles 
in which they should coincide.^ 

^Analogously, Type I theory on is T-dual to a Type IIA orientifold with D6-branes at angles [4]. 

''Notice that the two last refs. in [24] do not correspond to geometrical intersecting D-brane models, but 
to Type II orientifold constructions on Gepner points. The model-building techniques in these particular 
cases are, nevertheless, similar to those in [20]. For previous literature on Gepner model orientifolds see [26]. 

^Actually there is a small discrepancy in that we find that the F- functions appearing in eq.(77) of ref. [28] 



Due to the mentioned T-duality between the intersecting Dp-brane models on one 
side and magnetized toroidal Type I and Type IIB compactifications on the other side, 
all intersecting toroidal brane models constructed in ref. [4, 6, 18-24] admit a description 
in terms of the picture given in the present paper. In particular, the expressions here 
obtained may be used to compute the Yukawa couplings of the models in ref. [20] whose 
massless fermion spectrum is that of the non-SUSY SM or to those of the (local) = 1 
SUSY model in ref. [27]. With little changes they can also be applied to the computation 
of Yukawa couplings in the N = 1 orbifold models of refs. [25]. 

We would also like to emphasize that our results are also relevant for the phenomeno- 
logical extra dimension models considered in recent years [31]. In particular it has been 
suggested [32] that, if fermions of different generations have Gaussian- like profile in one 
extra dimension and are localized at distant points, one may obtain hierarchical Yukawa 
couplings from a one-dimensional overlap integral analogous to (1.1). Most of the specific 
models constructed [31,32] obtain chirality from solitonic scalar backgrounds. We find the 
approach in the present paper more effective since it gives rise to chirality, family replication 
and Gaussian-like profile in a natural and explicit way. 

The plan of the paper is as follows. In the next section we briefly discuss the ori- 
gin of Yukawa couplings in extra dimensional theories, leaving the details of the required 
dimensional reduction to Appendix A, and the details of SUSY compactifications to the 
Appendix B. In Section 3 we address the computation of the eigenstates of the Dirac and 
Laplace operators in T^" toroidal compactifications. These provide us with the wavefunc- 
tions of the light fields of the compactification, which are proportional to Riemann and 
Jacobi -iJ-functions and have a Gaussian profile. In Section 4 we generalize our results to 
the case with non-Abelian Wilson lines are present, in which the rank of the initial gauge 
group is in general lowered. Although quite interesting for rank reduction purposes, com- 
pactifications with non-Abelian Wilson line are technically much more involved and, up 
to some subtle points which allow to differentiate 9j from Ij/A, the final result for the 
Yukawa couplings is similar to the one obtained with the Abelian Wilson line case. The 
reader not interested in those details may safely skip this section. 

Armed with the wavefunctions of chiral fermions and light scalars we address in Section 
5 the computation of Yukawa couplings, by explicitly performing the integrals which mea- 
sures the overlap of three wavefunctions. In the case of SUSY compactifications we rewrite 
our results in terms of a superpotential and Kahler potential factors. This analysis is ex- 
tended in the following section to the case of D-branes of lower dimension (e.g. D5-branes) 
which are quite relevant for the construction of specific models. A detailed comparison 
with the results obtained in the T-dual models (involving intersecting D-branes) is given 
in Section 7, where the Yukawa couplings in both sides of the T-duality are seen to match. 
We also comment on the action that T-duality has on chiral fields, as can be deduced from 
this matching. As a simple application of the general results in this paper, in chapter 8 we 
present a 3-generation MSSM-like model and obtain the relevant Yukawa couplings using 
our formulas. Some final comments and discussions are left for Section 9. 



should have a 1/4 power instead of 1/2. This has also recently been independently pointed out in [30]. 



2. Yukawa couplings in Kaluza-Klein theories 



In this section we motivate the study of magnetized compactifications in order to achieve 
D = 4 chiral models from extra dimensions. We describe as well the general strategy that 
we follow to compute three-point functions in such models. We refer the reader to the 
appendix ^ or [1] for a more detailed discussion. 

2.1 D = W N=l Super Yang-Mills compactifications with magnetic fluxes 

Let us consider J\f = I D = 10 supersymmetric Yang-Mills theory, whose Lagrangian 
density is simply given by 

CsYM = -4^Tr {F^'^Fmn} + ^Tr {XT^'DmX} (2.1) 

where the trace is performed in the adjoint representation of a gauge group G, and M, = 
0, . . . , 9. The gauge group field strength Fmn and covariant derivative Dm are given by 

Fmn = OmAn - On Am - i[AM, An] (2.2) 
DmX = dMX-i[AM,X] (2.3) 

both the ten-dimensional vector Am and the spinor A transforming in the adjoint of G. 
Notice that the Yang-Mills coupling constant g has dimensions of (mass)^^ in D = 10. 

In order to obtain a D = 4 theory at low energies, we should consider the above 
theory compactified on a six-dimensional compact manifold Mq, so that we recover four- 
dimensional physics at energies below the compactification scale Mc- Indeed, in general 
the ten-dimensional fields Am and A admit a decomposition of the form 

A(x^ y-) = ^ xnixn ^ My"") (2.4) 

n 

AM{x^,yn = Y.^n,M{xn®^n,M{yn (2.5) 

n 

where n = 0, . . . ,3 and y™", m = 4, . . . , 9 stand for the non-compact and internal dimen- 
sions, respectively. The internal wavefunctions tpn, (pn M can be chosen to be eigenstates 
of the corresponding internal wave operator 

ipQtpn = rrintpn (2.6) 

A60n,Af = Ml^j^(j)n M (2.7) 

By applying the equations of motion, we find that the Dirac mass of the four-dimensional 

spinor Xn is given by m^, and so on. The lightest Kaluza-Klein replica then sets the scale 

of energy below which we recover a -D = 4 physics governed by massless fermions and 

light scalars. In practice, however, it turns out that plain compactifications on smooth A^g 

leave the gauge group G unbroken, the rank of G being usually too high to accommodate 

semirealistic interactions ^. Moreover, the computation of the Dirac index shows that such 

compactifications lead to non-chiral spectra in four dimensions [1]. 

^Whereas at the classical level any group G (or a direct product of these) is acceptable, at the quantum 
level anomaly cancellation imposes restrictions on the allowed gauge groups. 



Both unwanted features can be avoided by introducing non-trivial expectation values 
for the gauge field Am [7]. Indeed, since we are only interested in preserving Poincare 
invariance in the four non-compact dimensions, we are entitled to consider non-vanishing 
v.e.v.'s {Am{y)), m = 4, . . . , 9. On the one hand, the gauge group G will be reduced to a 
subgroup H C G commuting with the subgroup J which contains {Am{y))- On the other 
hand, a non-trivial gauge field modifies the Dirac operator and hence the computation of 
the Dirac index, and may introduce a chiral asymmetry that allows for a chiral massless 
spectrum [1,7]. We hence find that compactifications with non-trivial gauge fields {Am{y)), 
or equivalently, magnetized M.q compactifications with (Fmn) / 0, provide a natural way 
of achieving D = 4 chiral theories with reduced gauge group (See figure ||). 




Figure 1: Magnetized compactification. In the limit of large volume and diluted fluxes the field 
theory compactification with magnetic fluxes captures the underlying string theory physics. 

In addition, the introduction of a magnetic field in the compactification may not only 
lead to chiral matter but also to replication of chiral fermions, since the Dirac equation 
for the internal fermionic wavefunction pQip = may yield several independent degen- 
erate solutions, labeled by il^j{y)- In order to get canonical kinetic terms, these internal 
wavefunctions must satisfy 

I (fyi^jiy)^My) = Sjk (2.8) 

the same condition applying to bosonic wavefunctions. 

Finally, given the internal wavefunctions ipj, <j)k corresponding to the D = A chiral 
fermions and lightest scalars, it is possible to compute the Yukawa couplings between them, 
as an overlap between three wavefunctions. Indeed, the fermionic part of the D = 10 SYM 



action (2.1) contains a term of the form A - X - \, which upon dimensional reduction yields 
the Yukawa coupling^ 



^ijk 

IM 



[ i^f i;'^ c^l^ fate (2. 
Jm 



where fabc are the structure constants of the initial gauge group G. Notice that formula 
(|2.9|) provides us with the three-point function or normalized Yukawa coupling, i.e., not 



''See Appendix ^ for a more detailed discussion of the derivation of this formula. 



only contains the trilinear coupling of the superpotential Wijk, but also all the normaliza- 
tion factors coming from the Kahler potential/kinetic terms. Moreover, this expression is 
completely general, in the sense that we are not making any assumption on the holonomy of 
the compact manifold Mq or considering any particular embedding of the spin connection 
in the gauge group G. 

2.2 Models with fluxes on D = 6,8 dimensions 

The choice of D = 10 SYM naturally arises from considering the low energy effective action 
arising from heterotic and Type I theories, which are the simplest superstring theories 
involving gauge interactions. From the field theoretical point of view, however, we could 
consider e.g. the Lagrangian (|2.l| ) in D dimensions, D being even. The dimensional 
reduction scheme performed for D = 10 can be generalized for arbitrary D, obtaining a 
formula for the three-point function similar to (|2.9| ), but now with dim(7W) = D — A. 

It is in fact easy to build string based models where such D = 6,8 Lagrangians appear. 
Let us for instance consider type IIB string theory. Now, instead of compactifying it in 
a six-dimensional smooth manifold Mq, let us consider the possibility that presents 
a Zjv orbifold singularity, in such a way that in the vicinity of the singularity the metric 
can be written as T^" x C^~"/Z7v. We can then place a stack of N D(3 + 2n)-branes 
at such singular point, wrapping the T'^^ completely. In general, the worldvolume of 
D(3 + 2n)-branes in flat space will yield a massless sector containing aJ\f = lD = 4 + 2n 
U {N) gauge field theory, plus some extra matter transforming in the adjoint of this group. 
This extra matter is associated to the directions transverse to the D-branes, and signals 
the possibility of translating the D-brane in these directions. Now, since we are placing 
our D-branes at an orbifold point, these degrees of freedom are removed by the orbifold 
projection, and we are left with a much simpler massless spectrum, which still contains a 
£)-dimensional gauge theory. We can now turn on non-trivial magnetic fluxes F on the 
T'^^ which is wrapped by the D{3 + 2n)-branes. The final chiral fields will be given by 
the eigenstates of the Dirac equation on T'^^ which survive the Zjvf orbifold projection. 
Notice that the chiral fields live on the worldvolume of the D-brane and are hence trapped 
in 4 + 2n dimensions. Thus the overall wavefunction of the massless modes will have a 
non-trivial profile on T^" but will be a delta function in the rest of the dimensions. We 
then can apply the general formula (^j^) in order to obtain the Yukawa couplings between 
chiral fermions and scalars. 

Actually, this kind of construction is related by T-duality to certain intersecting DA 
and D5 brane models proposed in [6, 18], where the orbifold singularity was essential in 
order to get a D = 4 chiral spectrum. In fact, any of the semi-realistic models constructed 
in this intersecting D-brane picture can be translated to the magnetized extra dimension 
language, where the field theory techniques can be applied in order to compute quantities 
as, e.g., three-point functions. 

Let us illustrate these facts with a simple example. Namely, let us consider a geometry 
which is locally of the form x C/Zat, and place Na D5-branes wrapping and expanding 
the four non-compact Minkowski dimensions. The action of the orbifold group on the open 



Figure 2: General magnetized extra dimension scenario, involving compactifications with orbifold 
singularities. The gauge theory (D-branc) is just localized at the orbifold singularity. 



string spectrum is specified by a geometrical twist and the Chan-Paton action 

7^,a = diag (l;vo,e2-*l^i, . . . ,e2-^l^^^_i) , (2.10) 

where we have to impose X]£o^ -^a ~ If we choose a super symmetric orbifold twist 
then we obtain a D = 6 M = 1 supersymmetric spectrum, given by [6] 

Vector Multiplet n.=i UiK) 
Hypermultiplet ZZii^lK^') 

The addition of magnetic fluxes on will further break the gauge group to a product 
of unitary gauge groups. The techniques described in this paper can then be applied to this 
field theory. In particular, the Yukawa couplings are obtained by considering correlators 
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of light fields respecting the Z^v projection and performing the overlap integral over the 
2-torus. The final answer will then given by an expression like ( |1.2| ) for n = 1. 

Analogous D = 8 flux models may be obtained starting from Type IIB DT-branes 
wrapping a 4-torus and located at a Z^r singularity in the remaining two transverse dimen- 
sions. After addition of magnetic flux on the the gauge group is broken and chiral fields 
in bifundamental representations will appear. The massless modes will have non-trivial 
profiles on and again the Yukawa coupling will be obtained by computing an overlap 
integral over T'^ and imposing Z^v invariance. The final answer will again given by an 
expression like ( |l.2D but for n = 2. 

Again, note that these two classes of models are T-dual to the -D4 and D5 intersecting 
brane models of ref. [6] , and hence our formulae below provide us with the Yukawa couplings 
for these other classes of models. 

3. Toroidal wavefunctions I: Abelian Wilson lines 

In order to compute the Yukawa couplings from the formula (|2.9| ) we first need to have an 
explicit expression for the internal wavefunctions corresponding to massless modes in the 
low-dimensional theory. In the following we will compute such wavefunctions for toroidal 
magnetized compactifications. We will be particularly interested in how the wavefunctions 
depend on the Wilson lines involved in the compactification^ , since eventually the Yukawa 
couplings will have a rich dependence on those. 

We have decided to divide the computation of wavefunctions in two different sections. 
In the present one, we first consider the simplest possible case of wavefunction. That 
is a particle in a T^, charged under an Abelian gauge field, and in the presence of a 
magnetic flux. This example already illustrates the general form of the wavefunctions, 
which can be elegantly expressed in terms of Jacobi -iJ-functions. We then proceed to 
consider a more interesting example, which describes the wavefunctions in a magnetized 
compactification with gauge group breaking U{N) U{1)^. In this case chiral fermions 
charged under any couple of f/(l)'s appear. The results equally apply to a model with fluxes 
breaking U{N) — > Y\ U{Na) with ^ Na = N, in which case chiral fermions transforming as 
bifundamental representations appear. This is of course the phenomenologically interesting 
case in which e.g., a gauge group U{3) x U{2) x U{1) may be obtained, with chiral fermions 
transforming as bifundamentals. 

To be complete in the next section we will consider the more general case where non- 
Abelian Wilson lines are present. In this case the rank of the initial group will be reduced. 
Specifically, the final gauge group is of the general form U{Na) and J2a < ^- This 
more general case is also of great interest when considering a initial gauge group whose 
rank is much larger than six. It involves, however, a more technical computation and, 
for the specific purpose of computing Yukawa couplings enough insight can be gained by 
considering the 'Abelian' case. The reader not interested in the technicalities of the general 
case may then safely skip Section 4. 

*This Wilson line dependence has been usually neglected in the previous literature. It is, however, of 
central interest when considering Yukawa couplings. 



3.1 Eigenfunctions of the Dirac equation on 
3.1.1 Abelian gauge field 

Let us consider a flat two-dimensional torus ~ C/A, where A is a two dimensional 
lattice generated by ei = IttR and 62 = 2ttRt, r € C. The dual one-forms are defined as 
dx^{ej) = i,j = 1, 2. The metric is then given by 

ds"^ = gijdx'dx^ = Ihf.pdzf'dz^ (3.1) 

where dz = dxi + Tdx2. being a Riemann surface, any magnetic flux F solving Yang- 
Mills equations must be constant. In particular, let us consider an Abelian magnetic flux 
such that Jj.2 F = b, that is 

^ l;{dzA dz) = — ^ w, (3.3) 



ImT2 

where uj is the Kahler form derived from ( |3.2| ). Notice that, when expressed in terms 
of complex coordinates {z,z), no dependence of the area A = An'^R'^lm.T appears in the 
definition of F. This flux can be derived from the vector potential 

A{z) = -^lm{zdz), (3.4) 
2imr 

whose transformations under lattice translations are 

A{z + 1) = A{z) + ^Im dz ^ xi= 2ErTlm ^ 

(3.5) 

A{z + r) = A{z) + ^Imfdz ^ X2 = sETfl™^^ 

where we have deduced the corresponding gauge transformations A ^ A + dxi- We now 
consider a complex field (j){z) with charge q under the U{1) gauge field given by A. Its 
transformation under torus translations are given by (j){z) 1-^ exp {iqxi)<P{z), that is 

(f){z + 1) = 6*9X1(^)0(2) = exp 1^21^1™ z\ (l){z) 

, , (3.6) 

^{z + r) = e*'?X2W^(_2) = ea-p |i_^inirz} ())(z) 

Consistency of such transformations under a contractible loop in implies Dirac's charge 
quantization 

^ = M E Z. (3.7) 

We can now implement Wilson lines in this language. The simplest way is introducing 
a constant complex number C = Ci + TC2 such that 

«^^'"<-- + ^^' ,3.8) 

X2 = ^lmf{z + 



Notice that we can identify these gauge transformations as Wilson loops 



Xi = A{z), 71 (s) 



s 



X2 = §.y^^{z), 72(5) = ST 



s e [0, 1] (3.9) 



which implies that now we have a vector potential of the form 

A{z) = ^^Im ((z + Qdz) . (3.10) 
Im r ^ ' 

Any complex field (j) with charge q under such U{1), no matter in which Lorentz 
representation, will be described by a wavefunction that transforms in terms of X1jX2 C_| 
under lattice translations. We can express such wavefunction as 

0(^) = e^&-[(-+C)'l.e(z,z) (3.11) M 

where the function 6 must have the transformation properties 

e{z + i) = e{z) 

(3.12) 

Given the transformation properties of the function 9, we can expand it in a Fourier 
series along one of the real coordinates of z = x + ry. More explicitly, we can write 

0(x,y) = j;c„(y)e2™- (3.13) 

nGZ 

the second boundary condition in (|3.12| ) being translated into 

Cn{y + 1) = Cn+My) e-"^*^[^^^(l+2y)+ReC] (314) 

3.1.2 Dirac zero modes 

In order to compute the Dirac operator we need to specify a set of gamma matrices satis- 
fying the Clifford algebra. 

|pa^pfe| ^ 2,5'^^ (3.15) 
for we can choose the hermitian matrices 

In order to translate this algebra to the holomorphic coordinate frame, we proceed in two 
steps. We first consider the vielbein e 

9ij = eis'-dab, i.e., g = ■ 1 ■ e. (3.17) 

which in our case is 
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and which allows to express the Clifford algebra in the more geometrical coordinate frame 
by defining F* = e^F'*, being the inverse of the vielbein ( 3.17| ). These new matrices 
satisfy {F*, F-' } = 2g^^ , g^^ being the inverse of g in ( |3.17| ). In order to express this algebra 
in holomorphic coordinates we perform a further transformation / = {f,j^, f^) such that 

Ku = flfl9i, = fle<tfi^5ai„ i.e., h = (e ■ f)\e ■ f). (3.19) 

In our case 

r'=(]l] (3.20) 



1 r 



so that we finally obtain 



T-^i2.nr{liy r..(2,«)-.(°»). (3.21) M 



We also need to specify the covariant derivative 



Vi = di + LOi + Ai 



= di + ioif + AtT^ (3.22) 

where S*^' = i[F'=,F'] are the anti-hermitian generators of Lorentz transformations and Ta 
those of gauge transformations. Since we are in flat space we will take the spin connection 
cj^' to vanish. In the case of an Abelian gauge field on with Ta = —i we hence have 
= di — iqAi. We are finally able to compute the Dirac operator 

i^=|) +-^^ = i^F^Va + i5Zr"^'^ (3-23) 

a a 

which in our case is given by 



Let us now consider a two-dimensional spinor in 

^&(z,z-)= (3.25) 

transforming under U{1) gauge field ( |3.10| ) with charge q. Hence both wavefunctions 'ip± 
can be written in the form ( p. 11 ), ( p. 131) . Such spinor will contain a zero-mode of of the 
Dirac operator if it is annihilated by p, which implies 

D'ilj+ = ^ (d + ^^^lm {z + o] ■e{z,z) (3.26) 
\ Imr J 

= ^ (d- ^^^Im {z + C)^ • e{z, z) (3.27) 
V Imr / 
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Let us first consider eq.( |3.26| ). By using the Fourier decomposition ( 3.13| ) the zero-mode 
condition can be translated into the first order differential equation 



= -27rgM(Imry + ImC) + 2'irnT 



(3.28) 



kn being a constant. Inserting this solution into ( |3.13 ) we find that our wavefunction 

is given by 

^ ^^.gMiz+(V-^ . ^ kn e^™- (3.29) 
neZ 

thus being a holomorphic function up to a global function, as was to be expected from 
the general discussion in Appendix |^. Now, we can find the constant coefficients kn by 
imposing the recurrence relation ( ^.14 ) into ( 3.2^ ), obtaining 



(3.30) 



where Mn = Mn+qM are \M\ arbitrary constant coefficients. 

The existence of |M| independent coefficients signals the fact that the Dirac equation 
has \M\ independent solutions, each of them having a different wavefunction "0;^ j = 
0, . . . , l-M I — 1. Indeed, we can obtain such wavefunctions by splitting the summation index 
in ( |3.29| ) as n = r • M + j, and taking the overall factor J\fj out of the sum ( p.29| ). We then 
learn that we can express the final solution in a rather elegant way. Namely, as 



where is given by the Jacobi theta-function 



3 
qM 





{qM{z + C),qMT) 



(3.31) 



1? 



+l){u+b) 



(3.32) 



The fact that the final set of solutions ( 3.3l| ) satisfies the transformation properties 
defined by (|3.8[) can be easily checked by using the modular transformation properties of 
"(J-functions: 



1? 



(i/ + n,r) = e^™'^-?? 



{v + nr, r) = e-™'^-2™(^+^) . ^ 



(^,T) 



(3.33) 
(3.34) 



A parallel discussion can be done for ( 3.27 ). In order to summarize our results, let us 
define the function 



1_ 

N 



{Nu, Nt) 



(3.35) 
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Where the constants Nj will be soon specified. In terms of ( |3.35| ) we can define the 
wavefunctions 



(3.36) 
(3.37) 



where the star denotes complex conjugation. Is easy to see that 

D^^ = (g = +l) D\'ip{)* = {) (g = -l) 
Dt^^ = (g = +l) Z)(^i)* = (g = -l) 



(3.38) 
(3.39) 



We can interpret as the wavefunctions corresponding to left-handed fermions in 4D, 
while (V'+)* represent their anti-particles. Right-handed fermions would then be given by 
iIt'_ and their anti-particles by (ip'L)*- 

Notice the important fact that the solutions ( p.3§| ) and ( p.39| ) are mutually exclusive, in 
the sense that the theta function defining will only converge if M > 0, whereas ipi, is only 
well-defined when M < 0. Hence, by introducing a non-trivial flux M 7^ we automatically 
select one of the two chiralities of the two-dimensional spinor (|3l25| ). Moreover, we obtain 
several replicas of such chiral fermions, by means of the \M\ independent solutions of the 
Dirac equation. 



3.1.3 Normalization 

Once that we have found a basis of linearly independent wavefunctions, we proceed to 
express everything in terms of a orthonormal basis. This will allow us to have canonically 
normalized kinetics terms in four-dimensional reduced action. In terms of the internal 
wavefunctions we just found, this amounts to impose the following normalization condition 



dzdzijiiipi)* = 6jk, 



rp2 



(3.40) 



which is nothing but condition ( |2.8| ) for the particular two-dimensional case at hand. For 
the sake of concreteness, let us impose the normalization condition for ipj^. In terms of the 
wavefunctions ( p. 35 ) we have 



^{{^lY = V;i.'?M(r,z + C) • V'~'='-«^^(r,z + C) 
~ ^^■'«*^(t,z) • V"^'-5*^(f,z) 



qM 





{zqM, TqM) ■ ^ 



-k 

qM 





(3.41) 



{-zqM, -fqM) 





j 




k 






qM 


{z',t') ■ (^^ 


qM 















NjNk ■ e 



where in the second line we have get rid of the constant which will play no role when 
integrating over z. The last line is just the usual scalar product of theta functions, seen as 
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holomorphic sections of a line bundle over [33]. Indeed, integration over Kez imposes 
the condition j = k, and equality on the summation indices of the theta functions 



d{Rez) g-2vrM-[Imz]2/Imr . ^g27rQMImT(n+^)'^47rgAf(n+^)lm^ 

n 



It is now easy to integrate over Imz/Imr, since 



Im z 
Imr 



-27rgMImr(Ti+ 



j I Im 2 
qM 



^ flm.z\ -27rgMImTfn^ 



Imr 



j I Im 2 ' 
qM ^ Im T , 



(3.43) 



We thus find that, in order to satisfy (|3.40| ), that the wavefunctions ( 3.31 ) must be multi- 
plied by the normalization factor 



(3.44) 



where we have used the fact that, for ■0;^, qM = \M\. The computation for tp- give us the 
same result. 

3.2 Chiral matter eigenfunctions 

The previous section illustrates the general philosophy we will use in order to compute the 
internal fermionic wavefunctions in our magnetized compactifications. In the application to 
Magnetized Extra Dimensions, however, we need to consider a more general setup. Indeed, 
in general our effective field theory will be described by a U{N) higher dimensional theory 
broken to U{Ni), when turning on non-trivial magnetic fluxes. The chiral spectrum will 
arise from fermions transforming in bifundamental representations {Na,Nb)- Let us first 
consider that we have a two-torus with a magnetic fiux of the form 



Imr 



nib 



\ 



\ 



(3.45) 



where the rua a = a,b,c, . . . are different numbers. Again, by Dirac's quantization 
condition the magnetic quanta rua are given by integers, and it can be seen that turning 
on this kind of flux implies the gauge group breaking U{N) U{1)^.^ Let us then see 
how the previous computation of wavefunctions generalizes to this case. 



^See Section 4 for a more general, systematic discussion of these facts. 



3.2.1 Fermions in bifundamentals 



In order to compute the wavefunctions involved in a compactification with the flux 



( 3.45 ), it is enough to consider the simpler case 

(3.46) 



^ _ TTl I nia 



Imr \ nib I 

That is, a gauge group breaking C/(2) — > [/ x U{l)b- From the point of view of D-brane 
physics, this can be seen as two D-branes a and b, each with a magnetic flux in its internal 
world volume proportional to rua and nib, respectively. 

Dirac equation 

In the case at hand the gauge connection that can be chosen to be 

A, = ^Jr„ = -^f™"<^ + « (3.47) 

2Imr y mb{z + Cb) J 

Az = A^T^ = --^('^^^'+~^'^^ . (3.48) 
2ImT mb{z + Cb) J' ^ ' 

where = ^ * ^ ^ and Tb = . j are the anti-hermitian generators of C/(l)(j x 

U{l)b C U{2). The Dirac operator is again given by 

+ (3.49) 



The zero modes of ( ^.49 ) can be easily found. Let us consider a two-dimensional spinor in 



transforming in the adjoint of U{2). 

Again, ip^ will contain any zero-mode of the Dirac operator if it is annihilated by D, which 
implies 

= {nRr^ {di^+ + [A,, ^+]) (3.51) 
'd-&{z + Cab))c+ dD+ 



{ttR) 



-1 



where we have defined 



hb = ma - mb ^ 

Cab = {rriaCa " mbCb)/Iab 



(3.52) 
(3.53) 



On the other hand, the condition for tp- to contain zero modes is 

= {ttR)-^ (d^. + [A„ ?/;_]) (3.54) 

^ J /^^ , («-fi^(^-+«)B-\ „ 

\(«+5fe(^ + U))C- 80- ) 

which is consistent with the fact that we must consider ■0_ = ■0^, as ^' comes from a 
higher-dimensional gaugino. Notice that from ( 3.51| ) we deduce that A^,D^ have to be 
holomorphic functions, whereas C+ have to be of the form 

respectively, where ^{z) is an arbitrary holomorphic function, and is a z-independent ^ij 
normalization factor. 



M . ^±i^{z+uyim{z+u) . ^(^)^ (3.55) 



Gauge transformations 

In order to find the actual Dirac zero mode wavefunctions, however, we still have to 
impose them to be well defined in a gauge theory over T^, that is, to have the appropriate 
U (2) transformations under lattice translations. The transformations of an adjoint field 

A b\ . . 

^ ^ I are given by 

Eiz) ^ n..Eiz).nl=[^ ^ ^^^.(^)j.H(.).(^ ^ (3.56) 
where 

XfW = l^Im(z + Ca) 

a = a,b (3.57) 

We hence find that A and D are invariant under lattice translations, whereas B and C 
transform as 

B{z + 1)= e'^f^'^^Biz), B{z + t)= e'^2'(^)B{z) 

(3-58) 

Ciz + 1) = e-^^? (^)C(z), C{z + r) = e'^^^S (^)C7(z) 

where now 

Xf(z) = 5^1m(z + a) 

(3.59) 

Xt{z) = ^lmf{z + U) 
Notice that these boundary conditions are the same found in the previous subsection. 
Indeed, the transformation properties ( 3.58| ) can be understood in terms of the functions 
( |3.8| ), by making the substitution M lab and C, i— > Cafe- Finally, we must take q = 1 for 
the wavefunction B and q = —1 for C. We then find that the wavefunctions of ( |3.50| ) are 
given by 

B+ = ^,^'^^"{7, Z + Cafe), C+ = V^J'-^'^H^, Z + Cab) (3.60) 

and B± = (Cip)*. Notice that both wavefunctions in (|3.60 ) are again exclusive, and that 
B+ will then vanish unless lab > 0, whereas C+ will only be present for lab < 0. 
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On the other hand, we have found that A±, D± must be (anti)holomorphic and periodic 
under both lattice translations, for with the only possible solution is a constant function on 
T^. These constant wavefunctions are to be identified with the gauginos of the unbroken 
gauge group U (l)a x whereas the off diagonal entries of ^ are left and right-handed 

fermions transforming in the bifundamental representation of such gauge group. We then 
have 

i?+ < oo <^=^ lah > Left — handed fermions in (+1, —1) of U{\)a x ^7(1);, , . 
C+ < oo lab < Left - handed fermions in (-1,+!) of U{l)a x ^7(1)^ 

and thus we again obtain a chiral spectrum. 
3.3 Eigenfunctions of the Laplace equation 

The wavefunctions ( 3.60D turn out to be not only solutions the Dirac equation, but also 



eigenfunctions of the Laplace operator. The eigenvalues of such eigenfunctions, and hence 
the mass of the corresponding scalars, will in general depend on the Kahler moduli of 
compactification. In order to see this, let us compute the square of the Dirac operator 

(3.62) 

where the Laplace operator is given by 

2 

A = J^V^V^ = V? + Vi, (3.63) 

and Fzz by ( p.46| ). We hence get that the action of the Laplacian on such ab sector 
wavefunctions is 

same for (■0-^'='=^'''')*. 

These Laplace eigenfunctions turn out to be the ones with smallest eigenvalue. Ac- 
tually, the whole tower of eigenfunctions and eigenvalues can be recovered from them and 
the harmonic oscillator algebra [10,15]. Indeed, notice that for each ab sector we have 

A = iV + ^^, 
A 

N = D^D, (3.65) 

This reminds of the harmonic oscillator quantum algebra, which can be recovered by 
defining 



47r|I, 



ab I 



The eigenfunctions and eigenvalues of the Laplacian are then 



= (Dt)r^i,±/a6^ (3.67) 
= 27r^(2r + 1) > 0. (3.68) 



However, notice that this is not giving us the mass eigenvalues. From (A. 25) we know 
that the mass matrix is of the form 

M^=( ^ -^t^l (3.69) 



The eigenvalues of such matrix are 



A 



~Xr = Xr± iTT^ (3.70) 



which will give us one tachyonic scalar and the rest massive. The same tower of eigenstates 
will be present in the fermionic spectrum, the mass gap between the states being the 
same, but the lowest state corresponding to a massless fermion instead of a boson. This 
moduli-dependent spectrum describing a tower of massive fermions and scalars is the T- 
dual version of the 'gonions' described in [18], for the case of D-branes intersecting at one 
angle on T^. 

3.4 Theta functions as wavefunctions 

Let us pause our derivation of toroidal wavefunctions for a while, and try to gain some 
intuition from our results so far. In the previous subsections we have found that, after 
introducing a constant magnetic flux in a pure U{N) super Yang-Mills theory living in 
a two-torus, we obtain a broken gauge group G C U{N) and a chiral spectrum. The 
wavefunctions of such spectrum are given by either constant wavefunctions, which are to 
be associated with gauge bosons and gauginos of G, or by the non-trivial functions ip^{z), 
which will represent both chiral fermions and scalars in the dimensionally reduced theory. 
The mass splitting between fermions and bosons will be proportional to the flux density [3]. 

The chiral matter eigenfunctions are given by (|]3|) which, up to some normalization 
and exponential factors, are holomorphic Jacobi theta- functions. The exponential factor 
presents (up to phases) a gaussian behaviour on the torus coordinate y = S [0, 1]. On 
the other hand, the theta function dependence on z seems much harder to visualize. One 
can develop some intuition by plotting the square modulus of the wavefunction ip^'^ {t, z), 
associated to the 'probability density' p{z) = \il){z)\^ of finding a quantum particle with 
such wavefunction. 

We show such density in figure ^, for the particular case of a square torus (r = i) with 
one unit of magnetic flux = 1 and no Wilson line C = 0- We then find a square density 
which presents a gaussian-like behaviour in both axis of the T"^ x and y being, as well, a 
'periodic' function under lattice translations in both directions. The probability density of 
such function is peaked in (x,y) = (0,0), while it vanishes at {x^y) = (i, i). Of course, 
this is only in the particular case of C = 0, and these two points will be conveniently shifted 



Figure 3: Probability density p{z) = \ip{z)\'^ for a square with one unit of magnetic flux and 
vanishing Wilson line. The density p presents a symmetric, gaussian-like behaviour in the two axis 
of the X and y, being peaked at (0, 0). 

by varying the Wilson line. In general, the maximum and minimum of p will be placed at 
z = —Q and z = — Notice that the minima and maxima of a wavefunction density 
may be crucial in the final computation of a Yukawa coupling, since a 3-point function is 
given by a overlap of three such functions, and the minima or maxima of such can lead to 
enhanced or suppressed Yukawa couplings. 

Let us now consider a more generic case, namely when N ^ 1 and the spectrum of 
wavefunctions is composed of several replicas of the same chiral fermion/scalar. Let us 
choose = 3 (corresponding to 3 generations of the given fermion), which is moreover 
a phenomenologically interesting case, and plot the probability density = for 
J = — 1,0, 1. We show our results in figure g What we observe is that the three 
wavefunctions are similar, but shifted with respect to each other in the y = Imz/Imr 
direction, by units of 1/3 times the length of this radius. We find, moreover, that we have 
lost the symmetric gaussian-like behaviour. Indeed, each of the wavefunctions' density has 
a gaussian profile in the y axis of the torus, while it seems to be more or less constant in 
the X direction. Since we are considering C = 0) the gaussians are peaked at y = — j/S. 

At first sight it might seem quite striking that we have lost the x <^ y symmetry of 
figure ^. After all, there is nothing special in our problem regarding the y axis. A second 
thought reveals that this is just a matter of conventions. Namely, a matter of the choice 
of the particular basis which describes our space of wavefunctions. Indeed, a orthonormal 

^"Recall that in general j is an index defined mod \N\ so, for N — 3, j = —1 ^ j — 2, etc. 
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J = 

Figure 4: Probability densities p^{z) — \ip{zy-'^\'^ for a square with triplication of the chiral 
spectrum. The gaussian-like behaviour is now present in only one axis of the torus, namely in the 
coordinate y = Imz/Imr. The gaussians are similar to each other, and centered at the points 

y = -.1/3- 



basis for the vector space of wavefunctions with 'weight' N is giving by either [34] 



3_ 

N 





{Nv.Nt), j = l,...,iV (3.71) 



or 



{v,t/N), i = l,...,iV (3.72) 
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where N is given by ( |3.44| ) . These two bases are related by a discrete Fourier transform 



X 



k,N 



E 



(3.73) 
(3.74) 



Each of these two bases is suitable for describing the chiral spectrum of our compact- 
ification, being just related by a global unitary rotation in the space of wave functions. By 
plotting the densities of the elements of the alternative basis ( |3.72| ), we find a probability 
density where the roles of x and y have been interchanged, as figure shows. 





Figure 5: Probability density of an element of the alternative basis ( 3.72 ), p{z) = |x('2^)"''^P for 
a square T^. We have chosen N — 3, j = and C = 0, so that we are describing the same spectrum 
as in the figure ^. Notice that the squared density is a 7r/2 rotation of the wavefunction in figure ^. 



It turns out that the choice of basis il^^'^ or X'^'^ has a nice physical meaning from the 
point of view of string theory. We will unveil it when comparing our results with those 
found in Intersecting D-brane scenarios by means of T-duality. 

3.5 Generalization to T^" 

Let us now address how the previous computations generalize to magnetized compactifica- 
tions in higher-dimensional tori. We will first address the case where T^"' is a factorizable 
torus, and leave the more general case for the next section. Although the computations of 
wavefunctions becomes more technical, the final answer can be expressed as a product of 
n wavefunctions in T'^. 
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A factorizable 2n-dimensional torus can be decomposed as a product of two dimensional 
tori, that is 



r^" ~ X . . . X 



(3.75) 



In that case, both the vielbein e and holomorphic transformation / matrices are given by 
a direct sum of 2 x 2 matrices 



e = hnR^ 



(r) 



r=l 



1 RerW 
Imr(') 



^ ^ - ( 1 f{r) 



(3.76) 



r=l 



where r = 1, . . . , n labels each factor of in ( 3.75| ). As a consequence, the Clifford algebra 
in holomorphic coordinates can be reproduced by the following set of gamma matrices 



(2^i?('^)) 



1 
1 



10 0- 



(7^ (gl . . . (g) CT^ 



where o"^, o"^ are inserted in the r*'^ position, and we define 



a 



1 
-1 



2 




a 




2 



(3.77) 



(3.78) 



Now, let us consider the addition of a constant magnetic flux F in this background. 
We will consider F to be a (l,l)-form^^ in the complex structure defined by (3.75) and 
( |3.76D . We will again deal with the case where the gauge group breaking is given by 
U{N) U{1)^ . This implies having a magnetic flux of the form 



I rua 



TTl 



ImrW 



m 



(r) 



(r) 

rUc 



\ 



\ 



J 



(3.79) 



(r) 

where rria are n x N integer numbers, different for fixed r. 
B ifundament als 

In order to compute the wavefunctions of bifundamental fields, we again consider F to 
be a diagonal U{2) flux 



TTl 



m 



(r) 



ImrC'") 



m 



(r) 



(3.80) 



Following similar steps that the ones for the case and using the gamma matrices ( 3.77] ), 
we find that the Dirac operator is now given by 

ip= ij:,r^Dr-iEr'^''-Dl 



^^This condition is related to the hermitian Yang-Mills equations on gauge bundle compactifications. See 
Appendix H. 



Dr = a. + 



(3.81) 



2ImT('') 



The Dirac operator will act in a 2n-dimensional spinor, which can be expressed as a direct 
product of n 2-dimensional spinors. Each of the components of ^ is then given by V'ei,...,e") 
where e*^ = it. In this notation, the zero- mode equation is given by 



D^^^i^ =0 if e*^ = + 
dIiP.i ... = if = - 



Vr (3.82) 



On the other hand, the the boundary conditions for a field S transforming in the 
adjoint of U (2) is again given by ( |3.56| ), where now i = 1, . . . ,2n and 

X?.-i(^) = ^Im(z^ + C) 

^ a = a,b (3.83) 

X?.(^) = S^ImrM(z'-+a 
The wavefunctions of the 2n-dimensional fermions can be found by solving the differ- 



ential equations (|3.82| ) and the boundary conditions imposed by ( 3.56 ) and ( 3.83 ). As in 
the case, many of these solutions will be exclusive, and the ones which are non-zero will 
depend on the signs of the numbers I^J = m'^^ — m^\ The final answer is given by a tensor 
product of wavefunctions. In order to express it, let us first define the wavefunctions 

iC''^=^^'''^(r^'y + C^^ = (^^-^^''^^ (r^^\^^ + C))*. (3.84) 

where j^*) = 1, . . . , labels the 'Landau level' coming from each complex dimension, 

and ip^ ' ab is defined as in (|3.35D . The wavefunctions will then be given by 



/ const. n.'5.„ ..(,M) V'i 



ivn.5_^^,.(,(;))V'i const. 



(3.85) 



where s{I^^) = sign(/^2)- We hence find 2""! ^7(1)^ x U{l)b gauginos, associated to the 
constant wavefunctions, and lab = Ilr l-^afe I chiral fermions, associated to chiral fermions 
in the bifundamental representation. The wavefunctions of the latter are given by a simple 
product of wavefunctions of the form ( |3.6[l| ) and its complex conjugates. 

Laplace eigenvalues and masses 

Just as in the case, the wavefunctions ( p.85| ) not only represent zero-modes of the 
Dirac operator, but also eigenfunctions of the Laplacian. Is easy to see that the eigenvalues 
are now given by 

Again, one can recover the full spectrum of eigenfunctions and eigenvalues by consider- 
ing the superposition of n harmonic oscillator algebras, where the creation and annihilation 
operators Ur and at are defined in terms of Dj. and d}- as in ( 3.66| ). The spectrum of scalar 
particles is much richer than in the case, and in particular it does not have to be tachy- 
onic. We leave the derivation of the mass formulae issues for the general case treated in 
the next section. 



4. Toroidal wavefunctions II: non-Abelian Wilson lines 



In this section we consider the general case of introducing arbitrary fluxes, in particular 
those leading to a gauge reduction of the form U{N) Yii^iPi)^ with ^jPi < N. As 
we will see, this rank reduction occurs whenever we introduce non-Abelian Wilson lines in 
order to fulfill Dirac's quantization condition. In this case, there is a new technical issue 
when computing wavefunctions, which comes from the fact that a field transforming in the 
bifundamental representation {pa,Pb) is described by a x rif, matrix instead of a pa x Pb 
matrix, with > Pa- On the other hand, many of the details of the computation of the 
wavefunctions are similar to those in the previous sections, so we will be more sketchy in 
their derivation. — I 

4.1 Non-Abelian gauge groups ^j-j 

Let us first consider the case where we have a non-Abelian gauge group, say U{N). Both 
the magnetic flux F and the gauge potential A that we introduce transform in the adjoint 
representation of this gauge group. The Yang-Mills equations applied to the particular 
case of imply that any irreducible component of F must be proportional to the identity. 
Such irreducible magnetic flux is then given by 

m M , , 

Consider now a field $ in the fundamental representation of U{N). The main difference in 
this situation with respect to the Abelian case comes from the fact that the Wilson lines 
can be arbitrary elements of U{N). That is, the most general gauge transformation is of 
the form 

$(z) ^ ni{z)'^{z) = e'^^^^'^LJi'^iz) (4.2) 

where 

XiW = M^Im(^ + C)-liV, ^^^^ 

X2{z) = j^lmf{z + 0-lM, 

and LOi are constant elements of SU{N). 

Just as in the Abelian case, we can demand the gauge transformations to be consistent 
with the homology of T^, and to reduce to the identity when they correspond to a closed 
and contractible loop. That is, we impose 

n^\z + T)-n^\z + i + T)- n2iz + 1) • ni{z)<^>{z) = $(z) (4.4) 

Notice that this implies that the SU{N) part of ( |4.4| ) lies in the center of SU{N), that is 

• 1 ■U2-uJi = e^"^'^/^ -1^, k£Z. (4.5) 



After imposing this in ( |4.4| ), and by a similar argument as in the Abelian case, we find that 
M = k mod must be an integer, which is again Dirac's charge quantization condition. 
The LOi matrices can be chosen to be either [9, 13] 

oji = Q, uj2 = (4.6) 
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or 



M 



UJ2 = P 



(4.7) 



where we have defined 

/ 1 



P 



1 



1 



/l 



Q 



,2m /N 



,2ni{N-l)/N 



(4.8) 



are, of course, equivalent and describe the same physics at low 



Both choices (g^) and ( p 
energies. Each of them may be more convenient, however, for showing different aspects of 
magnetized compactifications. For instance, notice that the presence of the non-Abelian 
Wilson lines loi, uj2 imposes non-trivial constraints on the U{N) gauge potentials (or 
gauginos) transforming in the adjoint 



ioi ■ A. 



U]2- Afj,- 



(4.9) 



taking ( |4.6D and following [13], is easy to see that the constraints ( |4.9| ) impose A^ to be a 
diagonal N x N matrix, with only P = g.c.d.{N, M) independent elements. We thus see 
that introducing a magnetic flux reduces the rank of the gauge group as U{N) — > U{P). 

In fact, notice that if we have the particular case {N,M) = N ■ (1,5"), S" G Z, then 
we can take ui = uj2 = 1 in (|4.5| ), and the gauge group is going to be U{N), with no 
rank reduction. If we consider a D = 2 spinor ^ in the fundamental of U{N) there will 
appear S replicas of a chiral spinor transforming in the fundamental representation. The 
wavefunction of such chiral fermions will again be given hytp^'^ in (|3^ ). We then see that 
in the case where there is no rank reduction the computation of wavefunctions boils down 
to the ones performed in the last section, even if non-Abelian gauge groups are involved. 

Finally, in order to compute wavefunctions of chiral fermions and bosons transforming 
in bifundamental representations, as well as Yukawa couplings between them, the choice 



of non-Abelian Wilson lines (4.7) turns out to be more suitable. We will thus stick to that 
choice in the following. 

4.2 Fermions in bifundamentals 

Let us now consider F to be a direct sum of two such 'irreducible' representations. That is 



Imr 



rria -1 



nib-i 



(4.10) 



where + = A^. This can be seen as two magnetic fluxes Fa and Ff, over the same 
T^, both corresponding to two different gauge groups U{na) and U{ni,), and with magnetic 
quanta nia and mf, respectively. Actually, from the point of view of D-branes, each stack 
of D-branes corresponds to an 'irreducible' representation of F, so this system can be 
associated to two stacks of D-branes a and b wrapping T^, with multiplicities Ua and rib, 
respectively, and with a magnetic flux turned on each stack, each flux being proportional 
to nia and nib- 



As we have just seen, the low energy theory of such configuration will correspond to 
a gauge group U{pa) x U{pb), where pa = g.c.d. {ua, ma)- As noted above, if {na,rna) = 
na{l,Sa) then no non-Abelian Wilson lines appear and the gauge group will be given by 
U{na) X U{nh). The computations of the previous section will readily apply to this case, 
and the wavefunction of the chiral fermions transforming in (n^, n^) will be given by ( 3.6C| ), 



with the substitution lab ^ Sa — Sb- We then proceed to consider the cases where pa < Ua 
and some gauge reduction is involved. In order to simplify the discussion, we will suppose 
Pa = Pb = 1) the generalization to arbitrary rank being straightforward. 

Dirac equation 

Eq. ( |4.10| ) and condition (^^) now imply a gauge connection of the form 



2ImT\ ^-{2 + 4)1 



where Ta = {—i)lna and Tj, = (— are anti-hermitian generators of U{l)a x U{l)b C 
U{N). The Dirac operator is again given by ( p. 491 ). The zero modes of ^can be found by 
again considering a two-dimensional spinor of the form ( |3.5C1| ) but now transforming in the 
adjoint of U{N). The entry A± is no longer a number but a Ua y- Ua matrix, etc. We now 
impose the Dirac equation on the spinorial component which implies 



Di;+ = {TTRy' {di,+ + [A„ (4.13) 



, {d + + Cab)) B, 



= 



where we have now defined 



lab = -naTUb + mma ^ (4-14) 
Lb = lab/naUb (4.15) 
Cab = {nbrUaCa " nanibCb) / lab (4.16) 



which is a generalization of the definitions ( 3.52 ), ( 3.53| ) for > 1, and allows to dis- 



tinguish between the quantities lab and lab- This difference will turn out to be quite 
important. The integer number lab will again determine the multiplicity and chirality of 
the spectrum, and can be thought as a T-dual version of the intersection number from 
intersecting D-brane models^^. 

The conditions that the Dirac equation imposes on are also quite similar to the 
ones in the previous section, and can be obtained by making the substitution lab ^ lab 



and taking the general definition of Cab in (|3.54[) . Again we deduce that , have to 



^Notice, however, that we have defined lab with a relative minus sign respect to [6] 
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be constant matrices associated to U{l)a x U{l)h gauginos, whereas B+,C+ have to be of 
the form 

M ■ e=^*ST(^+C<'i')-i™(^+C<'6) .^(z), (4.17) 

respectively, where ^(z) is now an arbitrary holomorphic matrix-valued function, and M 
is a normalization factor. This time, instead of following the procedure of the previous 



section, we will consider the ansatze (4.17) and impose the U{N) gauge transformation 
properties that those fields must satisfy in order to be well-defined wavefunctions. It can 
be seen that both procedures lead us to the same final result. 

Gauge transformations 

The gauge transformations for a field ^{z) transforming in the {na,nh) representation 
of U{na) X U{nb) C U{N) are 



where 



$(z) ^ Q.1{z)^x){^\{z))^ = e^(>^^(^)->^^(^»w,'^$(z)(wj)t 
Xf(^) = ;i^lm(z + C.)-1„. 



X?(^) = ;^Imr(z + C„)-l. 



a = a,b 



^2 



P 



a = a,b 



(4.18) 



(4.19) 



(4.20) 



Hence, the boundary conditions for the components <pka,k^ of such bifundamental field are 
given by 



Hz+l) 



^{z + t) 



ei(x§-x|)p$(2)p-i 



(4.21) 
(4.22) 



. - ImT{2 + CQ^) 



where fee = 1, . . . , n^. Notice that ( 4.21 ) and ( [4.22| ) imply 



. J Im r(2 + Cfj^) 



(4.23) 
(4.24) 



Notice that these transformation properties will be satisfied by any field charged in 
the bifundamental representation (na, n^). Let us, however, focus on the zero modes of the 
Dirac equation. In particular, let us look for solutions of i?+ in ( |3.50| ), which must satisfy 



the ansatz ( 4.17 ) with a -|- sign. We find that the holomorphic matrix-valued function ^(z) 
is given by a theta function of the form 



(e^-(^))M=^ 



^ -I !_ 





{{z + Cab)Iab,Tlabnanb), 



(4.25) 
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where I = 1, . . . ,nani, has to be understood mod Ua, rif, respectively, and j = 1, . . . ,Iab- 
This solution is strictly valid and unambiguous only if g.c.d. (na, n;,) = 1, which is the case 
that we will consider in the following. We are thus able to express the bifundamental 
in terms of a linear combination of the following wavefunctions 

= J\f . e*i^(^+C«f')-i™(^+C«i') . ^^■'■^"'-(z), (4.26) 

and is easy to see that the hermitian conjugates expand a basis of wavefunctions 

for the bifundamental fields C_. 

Notice as well that ( 4.26| ) will only converge if lab > 0. In case lab < we will have 



fermions of opposite chirality, hence we should consider 



^j,4a and [^J'^bay (4 27) ["r] 

as the wavefunctions coming from fermions transforming in {na,nb) and {na,nb), respec- 
tively. 

Normalization 

The normalization condition in the case of bifundamental fields is given by 

J dzdz Tr (^<^''^-'^{^^^^-'')'f^ = 5ij, (4.28) 
On the other hand, notice that the boundary conditions (4.21) and ( 4.22| ) imply that 



+ 1) = (4-29) 

(4.31) 



As a result, if lab 7^ 0, we can compute ( 4.2^ ) by integrating 0^' "^^^(0;'; "'')* over a of 
complex structure naUbr}^ That is. 



/ dzdz TV = / dzd-zY,<p\MiiT 

dzdz<l^i0ii)* = {2luiT\Iab\nanby^'^ M^5ij (4.32) 



iT'^inariiT) 

from where we can extract the normalization factor J\f. 
Summary 

We have again found that the fermionic zero mode wavefunctions for chiral fields trans- 
forming in the bifundamental representation can be expressed in terms of theta functions. 
More precisely 

B+ = <^>^'^-», C+ = ^^-^-» (4.33) 



^^We are again considering g.c.d. (ria, nj) — 1. 
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where j = 1, . . . , \Iab\ and 



/ 2Imr|/, 



-1/4 



ab I 



I ^2 



(4.34) 



-fab "anb 





{{z + Cab)Iab,T'Iabnanb) 



These two solutions are exclusive in the sense that the theta- function series will not converge 
at the same time. Indeed^^, 



i?+ < oo <^=^ lab > Left — handed fermions in (+1, —1) of U{l)a x U{1) 
C+ < oo lab < Left — handed fermions in (—1, +1) of U{l)a x U{1) 



(4.35) 



The anti-particles of such wavefunctions are given by the hermitian conjugates of ( 4.33 ) 



4.3 Eigenfunctions of the Laplace equation 

Just as previously pointed out, the wavefunctions ( 4.34| ) will be also eigenfunctions of the 
Laplace operator ( ^.63 ). By a similar computation as the one in section 3.3 we obtain 



A A 



<j)i!±4b 



(4.36) 



same for ($-?'='=^''6)^ 

Notice that we recover the same eigenvalue than in the Abelian case of section 3, 
with the only replacement lab "-^ lab- We can carry out the quantum harmonic oscillator 
algebra and compute the eigenvalues of the mass matrix by making such substitution in 
the expressions ( 3.66| ) through ( |3.70| ). We finally obtain that the lightest scalar particle is 
given by a tachyon of mass 



"^tach 
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M 

A 



(4.37) 



This spectrum should match with the one obtained in the T-dual picture, at least in the 
limit of large volume A and diluted flux (small angles). By comparing both masses, the 
(approximate) analogue of the angle between two intersecting D-branes in the flux picture 
can be seen to be [16] 



T = - (*0"(>r9J - (an(i«t)) 

TT 



47r 



I. 



ab 



{A/a') 



(4.38) 



a quantity which only depends on the area of in string scale units and, in terms of the 
mathematical description of the magnetic flux as a bundle over T^, is given by the /x-slope 
of such bundle (see Appendix p| for the definition of //-slope) . 



^Recall that U(na) x U{nh) is broken by the flux to U{pa) x U{pb), where pa = g.c.d.{na,ma). The 



wavefunctions in (4.35) will in general be bifundamentals of such gauge group. In the particular case at 
hand we are considering pa = Pb ~ I- 
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4.4 Generalization to T^" 

Let us now address how the previous computations generalize to magnetized compactifica- 
tions in higher-dimensional tori. The most general constant magnetic flux associated to a 
U{N) gauge group is given by 

ft ' ' 

Fij = 2tt—^ ■ In (4.39) 

Naittj 

where riij = —riji. Here we represent T^" by the quotient R^"/A, where A = {x G R^^lx = 
nia^;n G Z^"}. We are then parametrizing the torus by the 2n-dimensional hypercube 
{x G R^"|0 < Xi < ai}, where ai = ||o'|| are the lengths of the 2n lattice vectors. 

The flux ( |4.39[ ) implies that the boundary conditions on a field <1> transforming on the 
fundamental of U{N) are of the form ] I ] 



<^{x) ^ ni{x)<^{x) = e*^^(^)a;i$(a;) (4.40) 

where 

and cJi are constant elements of SU{N). Consistency of the boundary conditions amounts 
to imposing 

nj^{x + ai) ■ n-^{x + ai + ttj) ■ Qjix + m) ■ J7j(x)$(x) = $(x). (4.42) 

This again implies that the SU{N) part of ( 4.42| ) lies in the center of SU{N), 

u;J^ ■ u;-^ ■ ijj ■ uJi = e^"^^-^/^ • l^v, Cij G Z (4.43) 

and that riij = Cij mod A^. Following [9,13], we consider SU{N) constant matrices P and 
Q such that PQ = QPe^'^'^/^ . The SU{N) part of the transition function can then be 
written as 

uJi = P'^Q'\ Si,tiGZ, (4.44) 
and the problem is reduced to finding Si, ti such that 

tiSj — tjSi = riij mod A^. (4.45) 

Such P and Q matrices can be taken to be (^^), the choices ( ^^ ) and ([4.7D being solutions 
of ( 4.45| ) in the particular case of T^. 

In a general T^" compactification, (4.45) will have a solution provided that all the 
higher Chern numbers are specified by the first Chern numbers of the fiux [16]. In the 
following we will assume that this is the case. Indeed, when dealing with constant magnetic 
fluxes, lack of satisfaction of (|4.45| ) means that the initial gauge group breaks as U (N) — > 
U{Pi) (instead of U{N) U{P)) after turning on F. The flux (B) can then be 
written as a direct sum of more 'fundamental' fluxes, each one giving rise to a gauge group 



in 
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U{Pi) [12]. For our purposes, then, we can just consider fundamental fluxes satisfying 
( |4.45[ ) and direct sums of these. 

In the following, we will generalize our previous results in order to compute wavefunc- 
tions of chiral matter fields in higher dimensional tori. We will follow the same kind of 
strategy as used for T^. We will first address the case where T^" is a factorizable torus. 
Although the computations of wavefunctions becomes more technical, the final answer can 
be expressed as a 'tensor product' of n wavefunctions in T^. We will then address the 
case of a general T^", showing that the wavefunctions can then be expressed in terms of 
Riemann theta functions. 

4.4.1 Factorizable tori 

Gauge group ^j-j 

Let us consider the factorizable background ( p.75| ) and the addition of a constant [-^ 
magnetic flux F on it. Without loss of generality, we will consider F to be fundamental 
in the sense described above and a (1, l)-form. This allows us to specify F in terms of 2n 
integer numbers (N^^\ M^^^), r = 1, . . . , n, such that [16] 



"-2r-l 2r = ^V^^^ • • • N^''-'^'^ M^''") N^''+^'l . . . A^W 



the rest of n^^i, vanishing. The components of the magnetic flux are then 

TTi 



(4.46) 



(4.47) 



ImrM A^M 

and the boundary conditions for a field <I> transforming in the fundamental are again given 
by Qr = exp(ixr) ■ ^r, where now 

X2r-l{z) = ^^^^lmizr+Cn-lN, ^^^^^ 

The action of the non-Abelian Wilson lines tOr on $ can be more easily described by again 
using a tensor product representation, now acting on the gauge group indices. Indeed, 
since A^ admits the decomposition (4.46), we can express the elements $ as 

k = l,...,N ^ fcM = l,...,iVM (4.49) 

we then have 

a;2T— 1 = ItvCI) ® . . . ® lAr(r-l) Q^'''' lj^(r + l) ... (8) l^in) 

W2r = lAr(i) . . . CS> lAr(r-i) -P lyvC-'+i) ® . • • ® 'i-^in) 



(4.50) 



where Q and P are the obvious generalization of ([4.8| ) to A^'-''-' x A^^*") matrices. After 
introducing such flux, the gauge group will be broken from U{N) to U{Y\j. P^"^^), where 
=g.c.d.(ArW,MW) [13]. 

^^In the mirror picture of intersecting D-branes, these fundamental constant fluxes correspond to D-branes 
wrapping T" submanifolds of T^" . 
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B ifundament als 



In order to compute the wavefunctions of bifundamental fields, we need again to con- 
sider F to be a direct sum of two fundamental fluxes: 



r-.T-^r — 

ImrM 



{'■) 



{'■) 



(4.51) 



where = Wr'f^a ■ The Dirac operator is now given by 

ip= fEfr^"^r-iE.™t 



/ ('■) 



\ 



(r) 

'2b ( _i_ /-r\ 



(4.52) 



and will act in the 2n-dimensional spinor, which can again be decomposed as the tensor 
product V'ei,...,e"- The zero more equation given by ( |3.82| ), we find that the solution can be 
expressed in terms of the wavefunctions 



-(r) Ar) ,-M / .f^. r(r)\t 



where j*-*^ = 1, . . . ,lj^fj, and ^^^"^^'-^ib is defined as in ( [4.34 ). The wavefunctions will then 
be given by 



(4.53) 



V'ei,...,e„ 



const. 



rJ 'I -"at I 



y 'I ao I 



Va.'^_,^,,(,M)(8).^t const, 
where s{I^^) = sign(/^^^). The tensor product in ( [4.54 ) is to be understood as 



'(8)*f''"-') 



n * 



^•J 'I ao 



(4.54) 



(4.55) 



Laplace eigenvalues and masses 

The eigenvalues of the Laplacian are in this case given by 



A (g)$i ' 



(4.56) 



and again, one can recover the full spectrum of eigenfunctions and eigenvalues by consid- 
ering the superposition of n harmonic oscillator algebras. We then recover the eigenvalues 
found in [15] 



(4.57) 



whereas the eigenvalues of the mass matrix are given by 



'{■">}, 



1/ 



Ml 



ab 



2s('^) + l)±2- 



fii) 



(2a'r' (j;i<),|(2.M + l)±20(^ 



(4.58) 



Note that the hghtest scalar excitations are obtained for s^"^^ = 0. In the case it is 
always tachyonic, reflecting the fact that = 1 SUSY configurations are not possible in 
this case. In the case the lightest scalar is either massless or tachyonic. Finally, in the 
T"^ xT^ X case the lightest scalar may be massive, massless or tachyonic, depending 

(r) 

on the values of the slopes 9app- In this latter case one recovers a SUSY spectrum if the 
lightest scalar is massless 

4.4.2 General tori 

Let us now consider the more general case where the 2n-dimensional torus is not necessarily 
factorizable. For simplicity, we will restrict ourselves to fields charged under Abelian gauge 
groups. That is, we set = 1 in ( [4. 391 ) for the rest of this section. 

A generic flat 2n-dimensional torus, T'^^ ~ C"/A, inherits a complex structure from 
the covering space C". Its geometry can hence be described in terms of a Kahler metric 
and complex structure as 

^ (4.59) 
dz^' = dxi" + rtfdy" 

where x^, y'^ E (0, 1), /Li = 1, . . . , n, parametrize the 2n vectors of the lattice A. The natural 
generalization of the Jacobi theta function ( |3.32| ) to this higher-dimensional tori is known 
as Riemann ??-function 



(i7|n) 



nieZ 



(rh+d) -CI- {m+d) g27ri(m,+a) ■ (u+b) 



(4.60) 



where a,b £ R", z G C"' and S7 is an n x n complex matrix. The transformation properties 
of such ??-function under lattice shifts are given by 



(i7 + n|n) 



2nid-n 



(z7|n) 



(z7+ n • n|ri) = e 



'TTin-d-n~2nin-{u+b) _ ^ 



{i7\Ct) 



(4.61) 



(4.62) 



where re € Z". These transformation properties are very suggestive. Indeed, inspired by 
( 3.35 ) we can construct the following wavefunction 

ij^'^{z, ft) =M- exp {i7r[N • • (N • Imn)-^ • Im [N • i]} • i9 



(N-z |N-n) 



^^See refs. [6, 16, 20] for a detailed discussion of the different possibilities in the T-dual language of 
intersecting D-branes. 



here AT is a normalization factor and fl. 



(N-zlN-fi) 



(4.63) 



Tnu- The transformation properties of this 



(4.64) 



wavefunction are given by 
provided that 

• N^^ e z 

• J • N G 

• (N-Imn)* = N-Imn 

and, of course, the series ( [4.63 ) will only converge if N • Imfi is positive definite. 

The natural candidate for the wavefunction of a field with charge (7 = +1 is then 
ip^'^(z + C) ^)) C now representing the Wilson lines. This wavefunction satisfies the differ- 
ential equations 

Da^P^'^{z + C^) = 0, Va 
Da = a^ + f f[N-(£+C)]-(Imn)-i' 



(4.65) 



and hence satisfies Dirac equation. If N • Im ft is positive definite, then it can be seen that 
these eigenfunctions satisfy the orthonormality condition 



by just fixing the normalization constant to 

Mn = (2"|detN| • det(Imn))^/^ • Vol(r2")-i/2 



(4.66) 



(4.67) 



where Vol(r2") = |det h| • det(Imn). 

In general, the integer- valued matrix N will encode the quanta of the magnetic flux. 
To see this more precisely, let us compare the transformation properties ( [4.64| ) with the 



1 2 4 

CO ljC cc cc 

ai ' 02 ' as ' 04 



transition functions (4.41) in the simple case of T^. By identifying 

' (x\y\x^y^) 

we recover the same transformations in both sides if we impose 

N = I I ^ 7^24 = ni3 = 



ni2 "-32 
ni4 77,34 



(4.68) 



(4.69) 



Now, as proven in [10], the degeneracy of states (i.e., the number of chiral fermions) 
is given by the absolute value of 

1 



Pf(n) 



e'^^^UijUki = det N 



(4.70) 
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at least in this particular case. In general, det N will give us the chiral spectrum obtained 
after turning on the flux: |detN| give us the degeneracy, whereas sign(detN) give us the 
chirality. 

Finally, notice that in order to have a well-defined wavefunction, the matrix N and ft 
must satisfy the following constraints 



• n24 = ni3 = 

• (N-Imfi)* = N-Imn 

• N • Im f2 > 



The first of this constraints is not such, since it can be satisfied by using the SL[A, Z) 
symmetry of T^. The other two can be understood in terms of supersymmetry, in particular 
from the requirement that F is a (1, l)-form (see Appendix]^. Indeed, in the case of T^", 
the sufficient and necessary conditions for (|4.39| ) to be a (1, l)-form are known as Riemann 
Conditions [33], which are 



-m* • Q • n 



H 

-H* 



(4.71) 



H being a n x n positive definite matrix. In our case, the matrices 11 and Q are given by 



n 



In n 
In 



n -n 

~ln In 



{lmn)-\ Q 



N* 



-N 



(4.72) 



and the Riemann Conditions amount to 



H 



(N • fi)* = N • n 

i N* -Imfi-i > 



(4.73) 



which clearly imply the constraints above. Finally, they also imply that, up to a phase, we 
can rewrite our wavefunction as 



N- (z + c)|N-n 



(4.74) 



5. Computing Yukawa couplings 

Once that we have derived both the fermionic and bosonic internal wavefunctions, and 
expressed them as an orthonormal basis, we are in position for computing the 3-point 
functions between them by using the general formula (|2.9| ). In this section we perform 
such computation for the toroidal compactifications previously considered. We will first 
focus on the simple case of T-^ and then generalize our results for higher-dimensional tori. 



5.1 Computing Yukawas on a 



In order to get non-trivial Yukawa couplings we need to start with three gauge factors, 
allowing for three different types of bifundamental matter fields. Let us compute the 
Yukawa couplings in the simplest case, namely magnetic flux compactifications in T^. In 
order to have non-trivial Yukawa couplings, we need to consider a flux of the form 



TTl 

Imr 



ria 



\ 



mt-. 



(5.1) 



mc-i , 



with Ua S JSS~^,ma G Zi, a = a,b,c. As explained in Section 4, the initial gauge group 
is broken to U{pa) x U{pb) x U{pc), where pa = g.c.d.{na,ma)- Notice that, with the 
definitions 



Ial3 = -nam/3 + npTTla, 
Ial3 = I a f^/ nan (3, 



(5.2) 
(5.3) 



the 'differences of fluxes' lap satisfy Iab + hc + Ica = 0. This implies that one {lapl is bigger 
than the other two. Let us suppose that this is the case for /^c, hence \Ihc\ = \Iab\ + \Ica\- 
This asymmetry will show up in the general formula for Yukawa couplings. 

We now have two possibilities, depending on the sign of Ibc- By the results of the 
previous sections, the fermionic wavefunction ^ is given by 




/ const. <I)*'^"i' 

const. 

y <j)i.^ca const. 

/ const. <I)J'^'"= 

<J)«,4a const, ^'^^^bc 

y const. 



if he < 



if he > 



(5.4) 



with ip- = (ijj^y, and where const. = A~^^'^ are the gaugino's wavefunctions. The chiral 
wavefunctions j = 0, . . . , |/q/3| — 1 have been computed in Section 4 and, in the 

particular case of Abelian Wilson lines, reduce to ip^'^°'f ■ M„^^„^, where ijj^'^'^f are the 
wavefunctions computed in Section 3 and Mn^^np is a Ua x n/3 matrix taking care of gauge 
quantum numbers. 



The general formulae ( |2.9| ), ( A.28| ) imply that the Yukawa couplings involving chiral 
massless fermions are computed by evaluating the integrals 



T2 



dzdz Tr {■(/;+ • [(/>_, '0+]} and 



T2 



dzdz Tr {ip^ ■ "0-]} 1 



(5.5) 



which are CPT conjugates of each other. Here (j)± are the wavefunctions of the bosonic 
fluctuations of the higher-dimensional gauge field Am, with helicity ±1 in the internal 
coordinates of T^. Notice that these are the only terms involving massless fermions allowed 
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by Lorentz invariance in the internal coordinates. By our previous results we find that (pj- ~ 
'iIj± and that it corresponds to the lightest (in fact, tachyonic) D = 4 scalar. Evaluating 
these expressions, we find that the Yukawa coupling involving left-handed fermions is given 

by 



ijk 



CTabc 9 Jt2 dzdz Tr <! '^''^'^ 



(5.6) 



where we have restored the dependence of the 3-point functions on the gauge coupling 



(Tabc 9 !t2 dzdz Tr {$*'^-' • 



constant g, by considering both (A. 28) and the normalization of the fields (A. 27). Here 
(^abc = ^^S^{labhclca) = sign(/aft/bc-^ca) IS & sign coming from fermionic statistics. From the 
point of view of = 4 physics, this term will yield a coupling between two chiral fermions 
of opposite chiralities, transforming in the {pa,Pb), iPa,Pc) bifundamental representations, 
and a complex tachyon in the {pb,Pc) representation.^^ 



Again, the computation of the integrals (5^) is technically simpler in the case where 
only Abelian Wilson lines are present in the compactification. The computation of the 3- 
point function is simplified by the use of -iJ-function identities. The case with non-Abelian 
Wilson lines has nevertheless the virtue that it can distinguish between two physically 
relevant quantities, which are the multiplicity of the spectrum, given by lap, and the slope 
of the flux (see Appendix]^), proportional to Iq,^. The differentiation of both quantities 
will turn out to be quite relevant when interpreting our results from the point of view of 
the effective action. 



5.1.1 Abelian Wilson lines 

Let us first consider the case which only involves Abelian Wilson lines. That is, we consider 
that (n^, m^) = nQ,(l, Sq) Sq, E Z. As shown in Section 4, turning on the flux ( |5.lD provokes 
the gauge group breaking U{N) — > U{na) x U{nb) x U{nc), Ua + rib + ric = and only 
Abelian Wilson line is needed for having a well-defined gauge connection. Moreover, the 



degeneracy of chiral fermions in bifundamentals {na,np) is not given by lap in (5.2) but 
rather by^^ Tap = Sa — sp. Finally, the matrix- valued wavefunctions <I>*'^°='3 reduce to a 



Ua X np matrix times the wavefunctions tf;^'^''!^ defined in (3.35). The Yukawas then read 



Yi,k = <Tabc 9 j^^ dzdz i^''^^^ {Z + U) ■ V'^'^^" {Z + Cca) ' (V''^'^^^ {z + Ccb)) ^ (5.7) 



where we have chosen I^c < for definiteness. The Abelian Wilson lines Cap are defined 
by ( |3.53|) , but with the substitutions nia ^ Sa and lap ^ ^ap- 



No more couplings are allowed from the choices made above and the action (A.l) in D = 6. In 
compactifications of higher dimensional theories, or those with a richer spectrum such as (2.11), however, 
more Yukawa between chiral fermions and scalars in different representations will appear. The computation 
of the 3-point functions in those cases will nevertheless be similar to the one we perform here. 

gives the degeneracy of chiral fermions before arranging them in bifundamental representations. lap 
is the topological invariant to be identified with the intersection number in the T-dual picture of intersecting 
D-branes. Both numbers agree when — 1, Vq. 
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Now, in order to compute the integral ( p.Tl) , we will make use of the 'addition formula' 
for theta functions, taken from [34], Proposition II. 6. 4. (p. 221), and which was crucial in 
the categorical mirror symmetry computations of [35]. This formula says 
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r 





N2 





{z2,tN2) = J2 ^ 
meZjv-^+jVj 



r+s+Nim 
N1+N2 





{zi+Z2,TiNi+N2)) (5. 



1} 



N2r-Nis+NiN2m 
NiN2iNi+N2) 





{Z1N2 - Z2Ni,tNiN2{Ni + N2)) 



This identity is particularly useful to our purposes, since the wavefunctions ip^'^ are 
proportional to -iJ-functions of this form. Indeed, 



^''^"H^ + Cab) ■ i^''^^^{z + Cca) = (2ImT)^/2 



,1/4 



(5.9) 



. gi^(2 + Ca6)Im(2 + Ca6)giTi^(2 + Cca)Im(2 + Cca)M 







3 





{(,Z -\- C,ca)'-^cai T^ca)j 



SO we can compute the product on the third line by using the formula ( |5.8[ ). We only need 
to identify 

r = i s = j 

Nl=Iab N2=lca (5.10) 

Zl = {z + Qab)1ab Z2 = {z + Cca)'Ica 

Let us first consider the case without any Wilson line, i.e., set Cafe = Cca = 0. We 
obtain 



1? 



lab 





{zlab,Tlab) ■ ^ 
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(z(-Jb,),r(-J6,))-t? 



'^^ab-^ca^bc 





(5.11) 

{0,T{-Iab1bc'^ca)) 



where we have made use of the fact that lab + ^ca = ^cb = —^bc- Formula ( 5.11 ) implies 
that 



^ab^ca 



'^bc 



1/4 



^i+j+^abm,Icb ^2) . -i^ 



Icgi-Igbj+Iablcam 
^ab^bc^ca 





(5.12) 

(0, r|JafcT{,c2^ca|) 



Now, the usefulness of these identities for performing the integration in eq.( |5.7D is that 
now the second theta function in the above expression no longer depends on the coordinate 
z and may be factored out from the integration. We are thus left we an integration over two 
theta-functions which may be easily computed by using orthonormality of wavefunctions, 
as we show below. Notice as well that ( ^.12 ) is a product of two wavefunctions, each one 
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with a 'weight' Tab and Z^a, expanded in a basis of a third class of wavefunctions, now with 
a 'weight' Icb = ^ca+^ab- Moreover this basis behaves under gauge transformations as the 
third wavefunction involved in the Yukawa coupling { f>.7\ ), more precisely as its hermitian 
conjugate. This kind of identity is totally general for magnetized compactifications, and 
comes from the simple fact that if we can understand a Yukawa coupling as an integral of 
three wavefunctions 

Y,,k= I Tr fC/3°</'L°'^U (5-13) 



M2 



with the trace performed over gauge and internal Lorentz indices, then the integrand 
must be invariant under both gauge and Lorentz transformations. In particular, (/)!^^ must 
compensate the gauge transformations of <Pap°<i^jj^- This is indeed the case, as can be easily 
seen, e.g., from the gauge transformation properties of {Na,Ni3) ■ {Nj3,N^) ~ {Ny,Na)^. 
This allows us to write the product of wavefunctions as^^ 



(5.14) 



where k runs over an orthonormal basis of wavefunctions transforming in the represen- 
tation {N^,Na)- The same facts apply to Lorentz indices. In the particular case of 
compactifications, we see from ( 5.12| ) that the coefficients are given by i?- functions. 
Now, by considering the orthonormality condition ( |3.4[lD , we can compute ( ^?7| ) to be 

1/4 
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(Tabc qA-^I'^ (2Imr) 
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'^^ab^bc^ca 
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^ca ^bc 





(0) T'l^ab'^bc^t 
/^ab 



{0,T\IabIbJca\) (5.15) 



So we find that Yukawa couplings are proportional to theta functions, as was already 
pointed out in the T-dual picture [27]. In order to compare both results, let us express the 
the ??-function characteristic in a more symmetric form. Notice that 
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'-be 
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^ab 



j-^bc ^ kl^d 
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Ibc 

^ca ^bc 



(5.16) 



where in the first equality we have made use oi i = k — j mod Tab, implicit in ( |5.15D , and 
in the last one we have made a redefinition of the indices j and k. This redefinition is 



^^As is stands the expression ( 5.14 ) is of field theoretical nature. Recall, however, that there is an 
underlying string theory in the whole construction, where the chiral fields 0' will be represented by vertex 



operators C. The expansion ( p.l4| ) is then understood as the field-theoretical version of the OPE 
in the underlying CFT (here uji is a world-sheet coordinate). 
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always possible if g.c.d. (lab, Ibc:^ca) = 1- In Section 7 we will see how this theta function 
characteristic matches with the result obtained in [27]. 

The inclusion of (Abelian) Wilson lines modifies the previous result to 



/2Imry/^ 

J-ijk — O-abc 9 I j 



'^bc 



1/4 



• e 



*T(2^abCabIm Cab+^bcCftcIm Cbc+^caCcalm Cca)/Im T 



• t9 



^ca ^bc 





but notice that, since Tab +^bc +Tca = 0, 

TabIca{C,ca Cab) — Qa^aXbc ~l~ Cc^cTab ~l~ Cb^bTcai 



(5.17) 



(5.18) 



Actually, it turns out that the whole Wilson line dependence of the Yukawa coupling is a 
function of the linear combination of Wilson lines ( ^.18| ). In order to see this, let us first 
express everything in terms of the following redefinition of complex Wilson line 



Co — SaCo) *^ — 6, C. 

The Yukawa couplings now read 



(5.19) 



where we have defined 
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C) T\-^ab-^bc-^a 
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k 



kjk -— + — + — 

-^ca -^6c 
Cafe = Ca - Cfe, etc. 

C ~ -^abCc ~l~ -^bcCa ~l" -^caCb 

Now notice that the exponential factor in ( 5. 20] ) can be rewritten as: 

gii- (CaSa/jIm ) /Im T 

where q, /? = a,b,c, and 5 is a symmetric matrix given by 

9 — l-Z^afe-^fec-^cal 



fee -^bc-^ca -^ab-^bc 



(5.20) 

(5.21) 

(5.22) 
(5.23) 

(5.24) 
(5.25) 



afe 



where we have again used the fact that Tab + 2^fec + ^ca = and that Xabl-bcX-ca < 0. This 
matrix is singular, having only one non-zero eigenvalue: 

1 ( ^^'^ 

9 ' 7= ^ca I — l-^ab-^bc^a 
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\Xab 



l^bA 

^ca I 

\Xab J 



(5.26) 



where k = +^bc ~^^ca- This impUes that the whole quantity in (|5.24 ) must depend on 
the combination C, above and, in particular, that we can rewrite ( ^.20 ) as 



where 
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H{C,t) = 27ri\lablbc^c, 
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(5.27) 



(5.28) 



We obtain similar results for Ij,c > 0. Note that the gauge coupling g in D = 6 has 
dimension of length, so that the Yukawa coupling in Z) = 4 is indeed dimensionless, as it 
should. 



5.1.2 Non-Abelian Wilson lines 



Let us now turn to the case where the flux in (5.1) satisfies pa = g.c.d.(nQ, m^) = 1 < Ua, 
so that non-Abelian Wilson lines have to be introduced and the total rank N = na + nb + Uc 
is reduced. The chiral fields are now expressed in terms of the matrix-valued wavefunctions 
^^Jai3 defined in ( [4.34 ), and the integral ( |5.6D is given by 
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dzdz 
/ dzdz 

JT^ 7_1 



where in the third line we have also assumed that g.c.d.(na, n;,, ric) = 1. 

The computation of this integral is harder than (^.6|), since we cannot use the theta- 
function identity ( |5.8D and the integral must be made by brute-force computation. Of some 
help is the fact that 



6ff='')*(z + r) 



^« + l,/ + l^i + l,/ + l*^^i + l,« + l^ ^'^1 



(5.30) 



And hence we can evaluate (5.30) by fixing I = and integrating over a torus of complex 
structure naUbncT, instead of performing the summation over /. We obtain the result 
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(5.31) 



where the Wilson lines have been redefined as = iTT-aCa, instead of ( 5.19 ). The exponen- 
tial prefactor is now given by 



H = 2-Ki\IabIbcIa 



Imr 



2-Ki ( laftCaftlm Cab + hcCbc^"^ Cbc + IcaCcal^ (ca ) /Im T 



(5.32) 



The expression ( 5.31 ) is quite similar to the one obtained in the case of Abehan Wilson 
lines. It is more general and contains more information, though, in the sense that we can 
now distinguish between the degeneracy of chiral states 7^/3 (intersection number) and the 
quantity I^p (/x-slope of the flux), which were both identified with Zap in the expression 



5.2 Higher dimensional tori 

The computation of the Yukawa couplings for 2n-dimensional tori can be easily deduced 
from the results obtained for T^. In the particular case that T^" has the factorizable 
geometry ( 3.75| ) the chiral matter wavefunctions are given by ( ^ ) or ( PI . That is, by 



j(r) r, 

a product of n wavefunctions of the form ^ ab or ($^ab )* obtained in T compactifications, 

(r) 

the choice of ^ or $* depending on the sign of I^^ . This implies that we can decompose the 
integral in as a product of n integrals of the form ( ^.61) . More precisely, each integral 



will be given by the analogue of ( |5.6D for the r^^ two-torus if 
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(r) 

and by its complex congujate if o"^^^ < 0. The final Yukawa coupling will be given by a 
product of n contributions, one for each T^, which are either of the form ( 5.31J ), either by 
its complex conjugate. 

We then find that the Yukawa couplings for factorizable T^" magnetized compactifi- 
cations are given by 
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whenever a 



where we must perform the substitution J*^*^ 1-^ J^*-* and C^*^ 
Here A^'^\ r^^^ are the area and complex structure of the r^^ T"^ in ( ^.75 ) and 

n 

Iab = \{lt^, 
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(5.36) 



is the total multiplicity of chiral fermions and scalars in the ab sector, etc., and the appro- 
priate labeling of them, with a different index i^^^ for each T^. On the other hand, 
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and \Ica \- Finally, 
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(5.40) 
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5.3 Yukawas in supersymmetric models 

Although the derivation of ( 5.33| ) is quite general, and in principle is valid for toroidal 
compactifications where super symmetry might be broken explicitly, it should be possible 
to understand it as a 3-point function inaAA=lZ) = 4 supersymmetric theory, at least 
in the cases where such construction can be achieved. The normalized Yukawa couplings 
that we have obtained should then fit in the general supergravity formula (see e.g.ref. [36]) 



(5.41) 



where Wijk is the corresponding trilinear coupling of the superpotential, K is the Kahler 
potential and Kab = dabdabK are the kinetic terms of the chiral fields in the ah sector, etc. 

There are indeed several examples in the literature of = 1 D = 4 chiral compactifi- 
cations realized as Type IIA intersecting D6-brane models [25], which are T-dual to Type 
I compactifications on a factorizable and with magnetic fluxes turned on.^'' Let us 
then consider the particular case of magnetized compactifications in T^, i.e., the case of 
n = 3, which involves D = \{) J\f = 1 super Yang-Mills compactifications with magnetic 
fluxes. This particular choice is not only relevant for Type I strings, but also for magnetized 
compactifications of heterotic strings and Type JIB involving D9-branes. 

In this case the coupling constant g is given hy g = e'^'^^l'^ ■ a'^/^, where 0io is the 
ten-dimensional dilaton and a' the string scale. The 3-point function then reads 
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and comparing this expression with (^HD we are led to the identifications 
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where in ( ^.44| ) we have neglected global phases and defined 
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Let us try reexpress the Kahler factors involved in the 3-point function in a more 
physical basis. Indeed, in terms of D = 4 supergravity fields, ( 5.44| ) reads 
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(5.46) 



^"Notice that, although aU of these supersymmetric models are based on orbifolds of which freeze some 
compactification moduli and impose some discrete symmetries to the open string sector, our computations 
and results are general and will equally well apply to these restricted geometries. 
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where the super gravity fields are defined as 



ReS = {2TT)-^e- 
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(5.48) 
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Notice that in ( ^.46 ) there is no dependence in ReT^^^ and the only explicit dependence 
on the Wilson lines comes from 

Y^l^^H^r)^ These Wilson lines are the vev's of the 
scalar fields in the adjoint of each gauge group, which belong to D = 4 = 1 chiral 
multiplets. Prom the point of view of D-brane physics, they are open string degrees of 
freedom, so their vev's compose part of the open string moduli space. Now, the Kahler 
potential for the closed string moduli has the well-known form 

3 3 

K = -log{S + S) - Zo(7 JJ +rW) - % JJ + (5.50) 

r=l r=l 

SO that we get an identity that the Kahler metric fields must obey in the flux side 
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6. D-branes of lower dimension 

As we have stressed in the introduction, J\f = 1 D = 10 SYM can be seen as the effective 
theory arising from type I and type IIB string theories at low energies. This fact has 
allowed to study D9-brane models with magnetic fluxes by using field theory techniques. 
From the point of view of type I or type IIB, however, D9-branes are not the only objects 
which yield gauge groups and/or chiral fermions in a general model. Indeed, both theories 
contain D5-branes and other objects of lower dimension than the D9-branes, which usually 
will carry additional gauge factors. Being of lower dimension than the D9-branes, these 
objects will appear as point-like in some directions of the internal compact manifold, while 
still expanding the four non-compact dimensions of the theory. These sectors of lower 
dimension are seen by the effective theory as small instantons, as is the case for heterotic 
strings. From the point of view of open string theories, however, these D-branes of lower 
dimension should be considered as natural as the D9-branes on which have based our 
previous analysis. In fact, it may turn out to be quite relevant for model building when 
trying to get a semi-realistic compactification, as the explicit model of Section 8 shows. 

Yielding a gauge group an chiral matter, it is fair to wonder which kind of Yukawa 
couplings do we obtain when considering D-brane models involving D5-branes. One would 
expect the general formula ( 5.42| ) to hold true, since under T-duality/Mirror symmetry 
both D9 and D5-branes are mapped to type IIA intersecting D6-branes, being on equal 



footing on this dual picture. Moreover, the Yukawa couplings between chiral fields in 
intersecting D6-brane models was computed in [27], and a general formula was obtained 
for every kind of D6-brane. 

In any case, it seems natural to wonder which wavefunction should be associated 
to an open string stretching between, e.g., a D9-brane and a D5-brane, specially in the 
dimensions where such D5-brane is pointlike. Intuitively, since the D5-brane looks like a 
localized source of RR charges/chiral anomalies in the D9-brane worldvolume, one would 
expect that the associated wavefunction is proportional to a delta function. In the present 
section we will try to further motivate this intuitive picture. 

In order to do this notice that, in the dimension where it is pointlike, a D5-brane 
can be formally thought as a D9-brane with infinite flux. From this point of view, the 
wavefunction associated to massless chiral fermion between a D9 with flux N^g and a D5- 
brane should be approximated by the wavefunction ( p. 60 ) in the limit lab ~^ oo, whereas 



that of an open string between a D9 and an anti-D5-brane should be recovered from the 
limit lab ~^ ~oo (here a labels the D5-brane and b the D9-brane). 

For the sake of concreteness, let us consider a D9-D5 system. We can then use take 
the formal limit ^ oo in ( 3.35| ). From the normalization procedure of wavefunctions, 



we know that 



im T J 



= A-h-'/^^e-^^h^^'-^y (6.1) 



where we have defined e ^ = 2ImrA^. Now, this last expression reminds of the family of 
gaussian functions which define a delta function 

d(x) = Lim,^o+^e-"4, (6.2) 

but with multiple images added (the sum over n) which is necessary to have a well defined 
function on instead of C. Indeed, if we take the limit — > oo, we recover 

n ^ 

This seems to suggest that in the N ^ oc limit, the wavefunction ( |3.31| ) should be propor- 
tional to (the square root of) d'j^2{Im{z + (i)/Imr). Notice, however, that for a D5-brane 
there is nothing special about the coordinate Imz on the whole problem^^. It thus seems 
sensible to consider the following wavefunction 

V^^^"" (r, z + C) ~ + (6.4) 



J.2 



^^As we have seen in Sections, this asymmetry may only appear when the hnear space of wavefunctions 
has dimension bigger than one, which is not the case. 



where some possible phase may be multiplying the whole thing. The former Wilson line C 
should be associated with the position of the D5 on the corresponding T^. 

In order to check that this is a sensible proposal, let us see which kind of Yukawa 
coupling we obtain by considering a system of two D9-branes, of fluxes nia and mb, lab = 
nia — mi) > and a D5-brane defined by the point The wavefunction coming from the ab 
sector will be the usual one, ip^'-^°-^{T, z + Cab)) whereas those coming from the 6c, ca sectors 
would yield a delta function, given by (^^). When integrating all of them we would then 
have 
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Notice that the whole expression can be put in the form ( 5.27] ), (|5.28 ), if we just define 

(6.6) 

which are indeed the intersection numbers arising in the T-dual construction^^. 



Cc =-Cc 
he = -1 
Ica = 1 



7. Comparison with intersecting brane computations 

We have emphasized several times along this paper the fact that magnetized compactifica- 
tions describe the same kind of physics as intersecting D-brane models do. In particular, 
D = 10 M = 1 SYM compactifications can be seen to capture the low energy physics of 
Type IIB string models involving D9-branes with magnetic fluxes. The latter are related by 
T-duality symmetry to Type IIA compactifications with intersecting D6-branes [37]. This 
pairing of theories can be extended to orbifolded/orientifolded versions of the above which 
involve, e.g.. Type I compactifications with fiuxes and orientifolded Type IIA models. 

Although describing the same physics, it should be clear that the techniques involved 
in the computation of = 4 physical quantities is quite different in both pictures. As a 
result, some elements of the D = 4 effective action may be much easier to compute in one 
side of the T-dual map than in the other side. This whole paper is meant to be an example 
of such important fact. Indeed, in [27] the Yukawa couplings of Type IIA intersecting 
D6-branes models were computed for toroidal compactifications, giving explicit formulae 
in terms of theta functions. The techniques used for such computation involved a sum 
over world-sheet instantons connecting intersection points. The Yukawas were obtained up 
to a global factor, h^u which encoded the contribution from quantum fluctuation of the 
world-sheet. The computation of such kinds of factors has been addressed in [28,29] where 
CFT vertex operator techniques were used. 

"^^ Recall that conventions on intersection numbers differ by a sign. Notice as well that in the particular 
class of models where D5-branes are involved the usual identity lab + he + Ica ~ does not longer hold. 
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In the present paper we have succeeded in computing the same quantity, the 3-point 
function, by using a different method. Indeed, the derivation of Yukawa couphngs in 
magnetized compactifications involve just field theory techniques, and no string theory is 
needed, at least in the appropriate limit of large volume and diluted fluxes. Notice as well 
that in this setup there is no 'classical' h^i and 'quantum' /iqu contributions to the 3-point 
function: everything comes from the overlap of three wavefunctions. 

In the present section we will check that the Yukawas computed from magnetized 
toroidal compactifications and D6-branes at angles match. More precisely, we will compare 
our present results with the ones obtained in [27] by using T-duality transformations, in 
order to extract the 'quantum' prefactor h^u and compare it with the one obtained in [28]. 
We will do such comparison in two ways, namely by performing a T-duality in one of the 
radii and then in the other one, obtaining the same result in both cases. As a byproduct of 
this computation we will gain some understanding of the action of T-duality on the chiral 
fields of the compactification. 

7.1 Intersecting brane Yukawas 

As was shown in [27], the computation of Yukawas in the context of Intersecting Brane 
Worlds in toroidal compactifications can be nicely expressed in terms of theta functions. 
More precisely, it was found that the formula for Yukawa couplings in a factorizable T^" 
was found to be 
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where the -iJ-function parameters are given by 
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and aabc = sign{lablbclea)- We are also supposing 
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In order to compare (7.1) with the results of Section 5, it is convenient to express this 
quantity in terms v.e.v.'s of complex fields. Indeed, notice that we can rewrite ( |7.1| ) as 
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with J^^^ defined as before and 
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In fact, expressions ( |7.6| ), ( |7.10| ) are vahd for ct^^|. = 1, Vr. In case we choose 



for any s we must substitute J^*^ i— > J^^'^ and fa 



abe 



-1 



, which amounts to taking the 



complex conjugate of the r*^ factor in ( |7.6| ). This is required by hermiticity of the effective 
Lagrangian, since changing the sign of the intersection numbers changes the chirahty of 
the fields at the intersection. Notice that we have found the same feature in the T-dual 
picture of magnetized D-branes. For simplicity, in the rest of this section we will consider 





Figure 6: Evaluating worldsheet instantons involving several boundaries will in general not 
yield a holomorphic function of the open string moduli Va = Oa + Jta, even in supersymmetric 
constructions, since their behaviour is asymmetric in 9a and ea- In the figure we illustrate the 
evaluation of the smallest instanton between three different boundaries. Notice that the area of the 
instanton will change, and so will the action, if we vary the position of the D-brane a, and we 
move it away from the point where all three intersect. The action will not change, however, if we 
stay at this point and arbitrarily vary the Wilson lines. See [27] for more details on these kind of 
compactifications. 

Finally, notice that ( [r.6| ) is not holomorphic in the open string moduli Va^. This may 
seem a little suspicious at first sight, since the 'classical' contribution is nothing but a 
sum ^ over holomorphic worlsheet instantons, and in heterotic models this provides 
a holomorphic function of the compactification moduli [38] . This seems also to be the case 
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for Type IIA compactifications on CY3 and with D6-branes wrapping Special Lagrangian 
cycles, at least if the open string instantons involve only one boundary [39,40]. When several 
boundaries are involved, however, we don't expect this to be the case, since the position 
and the Wilson lines which make up the complex field Va do not have a symmetric role 
when evaluating the world-sheet action (see figure |6|). Actually, evaluating holomorphic 
instantons (that is, surfaces calibrated by the Kahler form lo) only guarantees a expression 
holomorphic in the Kahler parameters, in this case J^^\ and not on the open string moduli 

(r) 

Va ■ Of course, in a supersymmetric compactification where we can understand Yukawa 
couplings as coming from a superpotential, the classical contribution h^i = (which 
encodes all the open string moduli dependence) must be holomorphic in all the fields up 
to a global term, that can be understood as a Kahler potential contribution to the 3-point 
function (5.41). The formula ( |7.6| ) is a good example of all these facts. 



7.2 Matching by T-duality 

In order to relate ( |7.6| ) to the results of Section 5 we need first to perform a T-duality 
transformation on the fields involved in the toroidal compactification. There are two dif- 
ferent ways of doing this, corresponding two the two different radii of T^. In each case the 
transformation properties are different, but we should obtain the same final result. We will 
perform both computations in the specific case of x x compactifications, but it is 
straightforward to generalize it to arbitrary number of tori. Before doing so, let us rewrite 
the 3-point function (5.42) in terms of more geometrical variables. First, instead of lap, 



let us consider the 'slopes' 
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which correspond to the spacing between particles in the tower of massive replicas ('go- 
nions') of chiral fermions and bosons. Actually, these quantities specify the harmonic 
oscillator algebra ( ^.66 ) present at each chiral sector of the theory. We will encode such 
slope dependence in the function 
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We can thus rewrite (5.42) as 
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and we have considered the four-dimensional dilaton (p^, defined as 
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J^^^ being the standard Kahler modulus of the r*'* two-torus^^ . 



7.2.1 Horizontal T-duality 

Let us first perform a simultaneous T-duality on the first radius of each T^, that is along the 
direction given by Im (tz) . We need to relate quantities by the T-duality transformation 



T ^ J 
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for each T"^ (for the sake of clarity, we will suppress the indices (r) from now on). We now 
make the substitutions ( [7.17 ) in (|7.6D , obtaining that the theta functions match exactly, 
whereas the exponential factors become 
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where ~ stands for equality up to a global phase upon exponentiation. We thus see that 
the part of the Yukawa coupling ( 7.13| ) depending on the open string moduli C precisely 
matches the classical worldsheet contribution to the Yukawa couplings computed in [27] . 
We are thus led to the identification 
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The idea is now to relate this prefactor ( 7.19 ) with the BCFT computation in [28]. Now, 
in [28] a quite similar expression was obtained as a prefactor. Namely, a square root of 
products of gamma functions given by 
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The quotient of Beta functions quickly decreases to one for small angles, as can be seen by 
using the definition of the Beta function and the Weierstrass form for the gamma function 
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which may also provide a way to measure the deviation of CFT computations from field 
theory results. We clearly recover the form of the prefactor ( [7. 11 ) for small angles. We do 
not, however, recover the square root of [28], but rather a 1/4 power instead. 

Despite the notation, ImJ^^^ should not be seen as the imaginary part of the complex field J^-^\ whose 
real part is not even a dynamical field in Type I models. Indeed, as seen in [22], it pairs with RR-field to 
form the complex scalar of an A/" = 1 chiral multiplet. 



Apart from the slope-dependent prefactor we are left with 
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The four-dimensional dilaton does not change under T-duality, so that we recover the 
prefactor 
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nicely matching the prefactors in [28]. 
7.2.2 Tilted T-duality 

We now consider the magnetized compactification again and perform a T-duality along 
the second radius of each T^, that is along the tilted direction given by Imz. After that, 
we get a intersecting D-brane model which has the same intersection numbers, and hence 
the same Yukawas, as the one we would obtain in the previous subsection. It will have, 
nevertheless, different wrapping numbers. 

Before doing that, let us rewrite the Yukawa couplings ( 7.13 ) in terms of the elements 
of the 'alternative' basis ( 3.72| ). 
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where Yijfc is given by ( 5.31| ). Notice that, if g.c.d.{Iab, ha Ica) = 1 then we can perform 
the following relabeling of indices 
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so that we have 
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The Yukawa couplings in the basis ( |3.72| ) are given by 
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Now we identify fields under the T-duality transformations 
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Making the substitutions ( |7.28| ) in ( [7.27| ), and using the transformation properties of the 
theta function 
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we get that 
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where we have defined 
We also get 
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Notice that the exponentials of ( [7.10D and ( [7. 321) match up to a phase. We are now left 
with 
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The identification of the slope-dependent prefactors works in the same way as in the pre- 
vious T-duality transformation, so after identifying them we recover 
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matching the previous result ( [7.23| ) after the transformation (|7.2g| ). 

Notice that both T-dualities lead us to the same 'quantum' prefactor in the intersecting 
D-brane setup. Moreover, by substituting the field theory quantity Q by the string theory 
analogue we recover a more symmetric expression in the angles. 
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where we have made the substitution 9^J = 1 — 0^^'' — OcJ implicit in ( [7. 201 ) 
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7.3 Chiral fields and T-duality 

The previous matching of Yukawa couphngs in magnetized D9-branes and intersecting 
D6-brane models suggest an intuitive picture of how T-duality acts on the chiral fields 
at intersection points. Indeed, notice that in order to match the Yukawas by means of a 
horizontal T-duality, we had to consider 3-point functions derived from the overlap of wave- 
functions of the form ( ^.71 ), whereas in order to compare results after a tilted T-duality we 



had to reexpress the Yukawas in terms of the wavefunctions ( 3.72 ). Recall that ( 3.71 ) and 
( 3.72 ) are the two canonical bases of wavefunctions on a T^. The matching above then sug- 



gest that there is a one-to-one correspondence between them and the chiral fields localized 



at D6-brane intersections after performing T-dualities (7.17) or ( 7.28| ), respectively. 

Let us reverse the point of view, and consider an intersecting D6-brane model. We 
can relate it to a type IIB model with magnetized D9-branes by performing either three T- 
dualities of the form (7.17) or those of the form ( |7.28| ). After the T-duality transformation 
the chiral fields, which previously associated to a pointlike intersection of D6-branes, will be 
represented by either the wavefunctions ( 3.711) (if we choose ( 7.17] )) or those in ( 3.72| ) (if we 



take ( 7.28 )). Notice that in both cases the wavefunctions of those chiral fields have their 
profile delocalized in the directions where the T-dualities have been performed), whereas 
their profile strongly depends on the other directions (see figures ^ and |5|). 
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Figure 7: Action of T-duality on chiral fields. T-duality maps the chiral fields at intersection 
points to wavefunctions defined on the whole compactification space. However, their probability 
density is delocalized in the directions of the T-duality and peaked in the transverse directions. 
Here we have considered the case of a horizontal T-duality acting on the wavefunctions in figure 4. 

This nicely matches with the intuitive picture that we have of T-duality of exchanging 
Dirichlet <-> Neumann boundary conditions, and hence changing the dimension of D-branes. 
Indeed, what in the type IIA picture was a D6-brane, localized in a 3-cycle of T^, after 
a T-duality becomes a type IIB D9-brane wrapping the whole compactification space^^. 
Actually, it is amusing to notice that these well-known facts about D-brane physics can also 
be understood in terms of delocalization of wavefunctions, this time those of the gauginos, 
which are given by a constant function localized in the worldvolume of a D-brane. 



Unless some dimension of tlie D6-brane is parallel to a T-duality direction, of course. 



7.4 The Kahler metrics of chiral fields revisited 



As we discussed in Section 5.3 in the magnetic flux picture, agreement with the supergravity 
formulae in the case of SUSY compactifications gives a constraint that the Kahler metrics of 
the chiral fields should obey. In the case of the intersecting D6-brane picture the analogous 
to eq.( 5.51 ) is 
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which yields 

3 



-3/2 



(2vr)3 r (e,)) r (etJ) r {e, 



r=l 



r 1 



n{r) 



r 1 



qir) 
^bc 



r 1 




(7.37) 



This seems to suggest a dependence of the chiral fields Kahler metric on the Wilson lines 
and twist angles of the form 
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where we have used a decomposition of H^'^^ analogous to the second line of ( ^.32 ). It would 
be interesting to compare these constraints with explicit expressions for the Kahler metrics 
of the chiral matter fields at intersections [30,41], as well as derivation of the Wilson line 
dependence of such metrics from direct BCFT computations, along the lines of [30]. 



8. A 3-generation MSSM-like Model 

As we mentioned above, the expressions here obtained may be used to compute the Yukawa 
couplings of the T-dual of the intersecting brane models in ref. [20] whose massless fermion 
spectrum is that of the non-SUSY SM or, in general, all of the toroidal models in [4, 6, 18, 
19,24,25]. For definiteness let us discuss here the T-dual of the model discussed in section 
4 of ref. [27]. This model has a spectrum quite close to that of a 3-generation MSSM and 
is simple enough so that we can display the explicit expressions for the Yukawa couplings 
without much complication. We leave a more systematic phenomenological analysis of the 
different models for future work. 

The model may be constructed as follows (see fig.|8|). Consider as starting point a U{A) 
D = 10 Yang Mills theory and let us compactify it on T'^ xT"^ x To obtain chirality 

^^This may be considered as a subset of a Type I string toroidal compactification. In fact, from Type I 
we should consider an SO{8) theory, broken to (7(4) by the addition of a non-trivial flux/ Wilson lines. We 
will not dwell in the details of the full Type I construction, though, since it would deviate us from the main 
point of this section. 




X 



T 



X 



T 



X 



T 



Figure 8: Structure of the 3-generation MSSM-like model described in the text. D9-branes with 
17(4) gauge group are wrapping the 6-torus and are subject to 3 units of magnetic flux in the second 
and third tori. One D5b{D5c) brane is wrapping the second(third) torus and is pointhke (black 
dots and crosses respectively) in the other two tori. 



we add twelve units of magnetic flux along the second and third torus (circles in the figure) , 
i.e., we have a flux of the form 

F,2,2 = j^U ^(ni'\mi')) = (4,12) = 4.(1,3) 
i^.3,3 = -j^l4 ^(ni'),mi=^)) = (4,-12) = 4 • (1, -3) 

As it was seen in Section 4, this flux does not break the 17(4) gauge group nor needs of 
the addition of non-Abelian Wilson lines. Hence, we can think of this configuration as a 
stack of 4 Z)9a-branes wrapping the 6-torus and each of them subject to 3 = 12/4 units of 
quantized flux on the second and third torus. We now add a couple of D5-branes: the first 
D5b is wrapping the second torus and is pointlike in the first and third tori (black dots 
in the figure). The second brane D5c is wrapping the third torus and is localized in the 
other two (crosses in the figure). It is well known that in Type I theory the gauge group 
in the worldvolume of n parallel 5-branes is USp{2n). In our case with two sets of isolated 
Db-hianes the overall gauge group will be [/(4)a x SU{2)b x SU{2)c, since USp{2) = SU{2). 
Now, solving Dirac equation on this background we will get bifundamental massless chiral 
fermions transforming under this gauge group as 3(4, 2, 1) + 3(4, 1, 2). We get three copies 
because of the three units of magnetic fluxes that we added. Thus in the end we get a 
Pati-Salam type of model with three generations of quarks and leptons. 

This is just a T-dual version of the model described in Section 4 of ref. [27] and more 
details about its structure may be found in that reference. Let us just mention that the 
U{1) in C/(4) is anomalous and becomes massive through a generalized Green-Schwarz 
mechanism in the standard way. Furthermore, one can break 5C/(4) SU{'3>) x U{\)b-l 
and SU{2)c U{1) by a judicious choice of Wilson lines on the worldvolume of the D5- 
branes, so that at the end of the day one is left with 3 generations of quarks and leptons 
and gauge group SU{3) x SU{2)lxU{1)y x U{1)b-l- The charged particle spectrum in 



this sector of the theory is supersymmetric if one chooses equal areas for the second and 
third tori Let us finally mention that there is a minimal Higgs sector like that in the 
MSSM if both L>5-branes sit on top of each other in the first torus (couple of nearby points 
in the first torus in the figure) All in all, the final spectrum in this subsector of the 
theory gets quite close to the form of the MSSM 

Let us now discuss the structure of Yukawa couplings in the above model. To simplify 
the formulae let us consider the initial Pati-Salam version with gauge group U{4) x SU (2) x 
SU(2), in which there is a single set of Yukawa couplings Yij{4.2, 1)*(4, 1, 2y{l, 2, 2). The 
only aspect which is a bit unfamiliar is the presence of D5-branes but we already showed 
in section 6 how to deal with them. Looking at fig. (|8|) we see that the first torus does 
not give us any fiavour structure since both 5-branes are pointlike and on top of each 
other (so that there is a massless Higgs field transforming like (1,2,2)). For the second 
and third tori we can apply the results in previous sections without any modification. In 
the present example we have (in the notation of Section 5.1) si' 



(r) 



ij^^ = a,P = a,h,c and T^p are the relevant numbers to be introduced in 



±3, so 



the general formula ( |5.27 ). They are given by 
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where we have also specified the relevant combinations of Wilson lines. In this case we only 
have two such combinations, corresponding to the U{1) inside ?7(4) and the second and 
third tori, with some contribution of the D5-brane positions Q and Q ■ The combinations 
are C^^^ = cP + CP and C^^^ = Ca^ + Cf \ and we can directly write down the Yukawa 
couplings in terms of them as 
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where e*^ is the D = W dilaton. One can easily show that these couplings are consistent 
with the ones obtained in [27] for a T-dual version of this model. There it was shown 



^''This is the T-dual of the SUSY condition in the case of intersecting branes, which requires equal complex 
structure in the second and third tori [27]. See Appendix^ 

^'^In other words, there is fj, mass term for the Higgs multiplets which is proportional to the distance 
between both _D5-branes in the first torus. 

^*It turns out that in the present model this subsector respects Af = 1 SUSY but the stringent conditions 
of RR-tadpole cancellation requires the presence of other branes which do not respect SUSY, so that the 
complete model is actually non-supersymmetric. Still its simplicity makes it a good choice to show an 
explicit result for Yukawa couplings. 

^^In order to relate both systems, one needs to perform a tilted T-duality described in Section 7.2.2 [19]. 



that in this simple model there is a single quark-lepton generation which becomes massive 
if the Higgs field gets a vev, which is a good starting point to reproduce the observed 
structure of quark and lepton masses. The existence of a single massive family is due 
to the factorization of the physics on the second and third tori of this particular model. 
Considering e.g., non-factorizable tori as discussed at the end of Section 4 would modify 
this. 

Note the following interesting point. Since the effective field theory from the heterotic 
50(32) string is also Z? = 10, = 1 SYM theory, a model similar to the above should be 
obtained starting from the heterotic string. The structure is identical, the only difference 
being that the SU {2)i, x SU{2)c gauge symmetries which came from Z)5-branes will come 
now from small instantons which are known to lead to simplectic groups. Thus the above 
model admits a (non-perturbative) heterotic construction involving small instantons. 

9. Final comments and conclusions 

In this paper we have computed the explicit form of Yukawa coupling constants in toroidal 
compactifications of I? = 10 SUSY Yang-Mills theories with constant magnetic fluxes. The 
results may also be applied to D = 6,8 extra dimensional models. The set-up studied 
is quite interesting, since it contains chiral fermions yet is simple enough so that one 
can explicitly obtain the wavefunctions of the light modes by explicitly solving Dirac and 
Laplace equations in the compact dimensions. This allow us to compute the Yukawa 
couplings as overlap integrals over the compact toroidal dimensions. Given the toroidal 
geometry is perhaps not surprising that both the wavefunctions and the Yukawa couplings 
obtained are given by products of Jacobi thcta-functions in the case of a factorizablc 
torus, and Riemann theta-functions in the general case. One would expect that in more 
complicated e.g. Calabi-Yau compactifications the Yukawa couplings will also be some type 
of automorphic forms. 

The class of models studied are T-dual to models of intersecting D-branes recently 
studied in the literature. Models with phenomenological interest have been constructed 
using that approach in recent years. We have found that, after appropriate redefinitions of 
the moduli and Wilson lines, the results obtained in both approaches agree in the dilute 
flux (small angle) approximation. This is interesting since both calculations are apparently 
quite different. In the flux case it is an exercise in Kaluza-Klein compactification whereas 
in the intersecting D-brane side is a stringy computation. 

The Yukawa couplings obtained depend on the complex structure and Wilson lines of 
the model. We have shown that the wavefunctions of chiral modes have a Gaussian profile 
in extra dimensions. Since the Wilson lines control the location of the maxima of these 
Gaussians, one can modify the results for Yukawa couplings by appropriately varying the 
Wilson line variables. Large Yukawa couplings should appear when the maxima of the 
three wavefunctions in the overlap integral coincide. On the other hand small Yukawa 
couplings should appear for wavefunctions with little overlap. This degree of freedom 
should be useful in order to reproduce the hierarchical structure of Yukawa couplings in a 
fully realistic model. Note that the Wilson line degree of freedom may be also understood 



as contributions from non-renormalizable Yukawa couplings. Indeed in the present class of 
models the Wilson lines correspond to the vacuum expectation values of singlet (or adjoint) 
complex scalars so that, expanding the Yukawa couplings on those fields one gets effective 
non-renormalizable couplings involving those scalars in addition to the fermions and the 
Higgs doublet. This kind of structure has been abundantly used in phenomenological 
analysis of Yukawa textures in the literature. 

As an example of the ideas studied in the paper we have briefly discussed a semire- 
alistic model with three quark- lepton generations. We leave for future research a more 
systematic study of explicit models which could perhaps be able to reproduce the observed 
quark/lepton masses and mixing. Note as a general remark that, due to the complex nature 
of both the tori complex strTicture as well as the Wilson lines, the Yukawa couplings are 
complex quantities and thus CP violation should be a generic property in a general model 
of this type. 
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A. Dimensional Reduction of N=l Super Yang-Mills 

Let us consider N=l Super Yang-Mills theory in D dimensions, D being even. Such theory 
is described by the action 

Sd = J d^w Cb + Cf (A.l) 

where 

Cb = -^Tr {F^'''Fmn] , Cf = ^Tr {It'^'DmX} (A.2) 

are the bosonic and fermionic part of the action, respectively, and M, N = 0, . . . , D — 1. 

The gauge group field strength Fmn and covariant derivative Dm are given, as usual, by (~ | 

Fmn = QmAn - On Am - i[AM, An] (A.3) 
DmX = dMX-i[AM,X] (A.4) 



where both the I?-dimensional vector Am and spinor A transform in the adjoint of the 
corresponding gauge group G. This action is invariant under supersymmetry as well as the 
gauge transformations 

Am Am + SmO + i[e, Am] (A.5) 
X^X + i[9,X] (A.6) 

6 being an arbitrary function of w taking values on adjoint of G. 

For the sake of concreteness let us choose the gauge group G = U (N). The Lie algebra 
basis of such group can be chosen to be (f^a)} = ^ai^aj, {e-ab)ij = ^ai^bj- We can then 
expand the fields in the adjoint Am and A in terms of such basis of generators as 

Am = Bm + Wm = BljUa + W^Cab (A.7) 
A = X + ^ = x'^Ua + ^"'e„b (A.8) 

Hermiticity of the U{N) generators imposes Bm and x to be real and W^^ = (Wlf)*, 



\^ab _ ^i^bay ^ Substituting ( |A. into the fermionic Lagrangian £f we find 

Cf = cP + Cy' + Cf' (A.9) 

where 

cP = -\tv {^r^dM^-i^r^'[BM,^]} , 

JOy' = {^T'^'iWM,^]} , (A.IO) 

Cf' = ^TV {xT'''dMX-ixr^'[WM,^-i^T^[WM,x]}, 
and performing an analogous computation with the expansion ( |A.7D and £b we get 

Cb = cf + cf + Cb' (A.ll) 
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with 



(2)' 



B 



(4)' 



1 



Tr {DmWnD^^W^ - DmWnD^W^^ - iGMN[W^,W^]} 



B 

Cb' 



Ag 

i 



■Tr {[Wm.Wn][W''\W'']] 



(A.12) 



1 



^^^Tr {{DmWn - DnWm)[W^ .W]] - ^Tr [GmnG^""] 



where we have defined 



(A.13) 
(A.14) 



Gmn = QmBn — QnBm 
DmWn = BmWn - i[BM, Wn]. 



Now, let us expand all fields in terms of their Lie algebra components as in ( [A.7| ), 
dAj). We find 



(2) 
B 



1 

1 

2? 



(A.15) 



c 



(2) 



F 



2r 



2^2 



qiabYM-^yca^bc 



where Dm'^ is defined as in ( A.14| ). 



Notice that, up to now, we have done nothing but reexpressing equation ( A.l ) in terms 



of the new fields B1^{w), W^f{w), x"(^) and ^ (w). Let us now compactify this theory 
in a — 4 dimensional manifold Md-a performing the dimensional reduction from D 
to 4 dimensions in two steps. First, let us decompose the D-dimensional fields B'^{w) 



and W\ 



ab. 



N 



W 



into the usual components B^{w), Bf{w), (w) and (w) n = 0, ...,3, 
i = 4,...,D — 1. Since we are interested in maintaining Poincare invariance in the fx 
coordinates, we are free to give non-vanishing vevs for Bf^{w) and W^^{w): 



Bt{w) = {Bt){y) + Ct{w) 
Wf{w) = {Wf){y) + ^t{w) 



(A.16) 
(A.17) 



Notice that these vevs can only depend on the compact coordinates y if we want to pre- 
serve four-dimensional Poincare invariance. These non- vanishing vevs will generically break 
gauge invariance. These expectation values correspond to the turning of magnetic field in 
the compact dimensions. In the following we will consider {Wf'){y) = 0, that is, we 
will restrict ourselves to Abelian fields in the compact space. Moreover, the case we will 
ultimately be interested of is constant magnetic fields. 
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Since our interest in this paper is to compute Yukawa couplings in the lower dimensional 
theory, we focus on the terms in ( A.15| ) involving D = A fermions and scalars. Such scalars 
are given by the fields Cf{w) and <^f'{w), so let us rewrite everything as 

.(2)' 



c 



B 



A2) , r(2) 

'-B "T '-B 



2g \ / \ 

2g L 



iab\ 



+ c 



,(4)' 



A^) , A'i) 

'-'B "T '-B 



1 

2^ 



^ab^hc^icdf^jda _ ^ab ^bc cd ^ida 



Cy' 



29 

Cy + Cy 



1 

2^ 



ab 



+ A 



(A.18) 



where the Cs contain terms irrelevant for our subsequent discussion and we have defined 
the 'average' covariant derivative by Di = di — ig{Bi). Notice that '^"'^ transform in 
the bifundamental representation of the (D-dimensional) gauge group U{\)a x ?7(l)b, hence 
this derivative acts as 



Di^f = di^f - i{B^)Wf + i{Bl)W° 

,(2) 



(A.19) 



same for and contain all the possible quadratic terms^'' in $j and that 

can give rise to effective mass operators for the dimensional reduction of these fields. 

The second step is to expand the D dimensional fields on a basis of eigenstates of the 
corresponding internal wave operator: 



n 



(A.20) 
(A.21) 



and so on for the rest of the fields, satisfying 



ab 



(A.22) 
(A.23) 



where r(4) = iV^T'^T^T'^. We are choosing the internal wavefunctions il)'^,(f)'^^ to be di- 
mensionless. The eigenvalues • are directly related to the four-dimensional mass 

of the fields Xn'') V'n^r Indeed, by applying the equations of motion we find 

(A.24) 
(A.25) 



-mnxf 



2^| {G^ 
IM 



°In the case of non-zero vevs for non-Abelian internal gauge fields < W >iab (y) 7^ more terms like 



e.g. DmWn[W^'',W^] should be included in c'^'' and £ 



,(2) 
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Notice, as well, that in order to get canonical kinetic terms the {D — 4)-dimensional 
fields must satisfy [2] 



g-' / d^-^y = ^acS^d (A.26) 

Jm 

g-^ [ d^-\ r^ivi'r^iy) = SacSbd (A.27) 
Jm 

Functional integration of the massive fields would yield a four dimensional theory 
including just the lightest modes. The massless four dimensional fields would be given by 
U{na) gauge bosons A'^, the gauginos A^, the scalars in the adjoint cf, and bifundamental 
spinors Xo^- There will also be bifundamental scalars (pf^, which may be massive, massless 
or tachyonic depending on the details of the compactification. In a semi-realistic scenario, 
the D = 4 gauge group obtained after dimensional reduction is identified with Standard 
Model gauge group, the bifundamental fermions with the SM quarks and leptons and 
one bifundamental scalar with the Higgs boson. Finally, the several replicas of each field 
that may exist due to multiplicity of zero modes are identified with the different particle 
generations. 

One of the main motivations of this paper is to compute the generic form of the Yukawa 
couplings for the lightest D = 4 fields in the class of compactifications described above. 
Using ( A.igj ) and ( A. 21 ) we find that such couplings are given by 



- I d'x xf V>ji • 1^ d^-'y ct>7, r V^^^l (A.28) 
where the {/, J, index the replicas that may exist for each of these fields. 



B. Fluxes and supersymmetry 

The above discussion is also general in the sense that it does not depend whether the 
D = 4 low energy-theory is a supersymmetric field theory or not. When coupling SYM 
theory to gravity, however, it may be useful to consider compactifications where N = 1 
supersymmetry is preserved in the effective theory, at least at the perturbative level. The 
amount of supersymmetry preserved at low energies depends on the geometrical details of 
the compactification. In the following, we briefly review the conditions for supersymmetric 
magnetized compactifications. We refer the reader to Chapter 15 of [1] or ref. [42,43] for 
more detailed discussions. 

B.l Hermitian Yang-Mills Equations 

Consider D = 10 SYM theory compactified in a 2n-dimensional manifold coupled 
to = 1 supergravity. Under the assumptions of H = d(j) = 0, the conditions for a local 
unbroken supersymmetry amount to having a covariantly constant spinor rj. This in turn 
implies that the compact manifold A^2n is a Ricci-fiat Kahler manifold, i.e., a Calabi-Yau 



n-fold CYn. On the other hand, a non-trivial Yang-Mills gauge field A must have a field 
strength F = dA which obeys 

6\ = r^FijT] = (B.l) 

where A stands for a ten-dimensional gaugino. If A^e is a complex manifold (as the existence 
of the covariantly constant r] would require) then we can rewrite (|B.1|) as [1] 



Fab = F-^-, = (B.2) 
a^'F^i = (B.3) 

where g^^i is the hermitian metric on A^2n- These two conditions are quite strong. For our 
purposes we may consider a slightly less constraining system of equations, which comes 
from coupling D = W SYM to = 2 supergravity. Eq. ( |B.1| ) then generalizes to 



SX = r^FijT] + 77' = (B.4) 

where rj' is the other covariantly constant spinor coming from extended supersymmetry. 
Eqs. (Q, then relax to 

Fab = = (B.5) 
g-'F^j, = c-l (B.6) 

where c is a constant which encodes which M = 1 subalgebra of = 2 is preserved by 
the SYM theory. Eqs. ( [B.5[ ) and (|B.6| ) are known as hermitian Yang-Mills equations, and 
imply that the gauge field A solves the super Yang-Mills equations, hence giving rise to 
a supersymmetric gauge theory upon dimensional reduction. If TWg is a Kahler manifold 
(which is again implied by the existence of a covariantly constant spinor), we can rewrite 
( |B:6D as 

F A = c • 1 (B.7) 

where u) is the Kahler form of Mq. 

Equation ( [B.5| ), provides a generalization of a holomorphic function on A^2n- Indeed, 
notice that ( [B.5D implies F^i = i[Da, -D5] = i[da — iAa, — iA^] =0, so if we assume that 
the gauge field A is hermitian, we can write it as 

A-a = -id-aV ■ V-\ Aa = -i{daV^-^) ■ V\ (B.8) 

where y is a matrix-valued function on the coordinates (2", z"). Then, if we define a field 
/ to be "holomorphic" by satisfying 

D-af = (B.9) 

then is easy to see that such field is of the form / = Vg^ where g' is a holomorphic function in 
the usual sense, i.e., dag = 0. In geometrical terms, the gauge field A can be understood as 



a vector bundle E over M2n- Now, a vector bundle satisfying (B^) is a holomorphic vector 



bundle, in the sense that their transition functions can be chosen to be holomorphic [1]. 



The converse is also true. Hence, considering Yang-Mills fields satisfying ( |B.5D amounts to 
studying holomorpliic gauge bundles. 

On the other hand, eq. ( |B.6| ) is related to the notion of /x-stability by the theorems 
of Donaldson and Uhlenbeck-Yau. Indeed, let us consider an irreducible complex vector 
bundle E over a compact Kahler manifold M2n, and whose field strength is given by 
F = dA. The degree of E is defined by 

deg{E)=f ci(S) Aw""^ = / trFAw""^ (B.IO) 

where ci{E) is the first Chern class of E. The relevant quantity is the slope of E which 
is defined as fJ.{E) = deg{E)/rank(E).^^ Notice that fi depends on the Kahler moduli of 
M.2n, being independent of the complex structure moduli. A bundle E is said to be stable 
if for any subbundle S C E we have /x(S') < fJ-{E). What the theorems above show is that 
any bundle E satisfying (|B.5|) (that is, a holomorphic bundle) and being ju-stable satisfies 



( B.6 ), and the other way round. 
B.2 D-brane interpretation 

In general, we will consider a magnetic flux which is of the form 



F = \ Fb I (B.ll) 



Each box Fa, a = a,b, . . . represents an irreducible component of the vector bundle (that 
is, a subbundle). Turning on the magnetic flux F will break the initial gauge group G to 
a direct product of smaller gauge groups Gi C G. Generically, each of these smaller 
gauge groups Gi will be an Abelian U{l)a factor associated to the component Fa- If we 
consider a compactification where Fai = Fa, i = 1, . . . , Na then we may have the gauge 
enhancement Yia^ ^ U{Na).^'^ 

This picture reminds of the gauge theory dynamics associated to D-brane physics and, 
in fact, it turns out that D-branes provide a nice physical realization of the mathematical 
results previously stated. In order to see this, consider type IIB string theory compact- 
ified on A4q with D9-branes filling up the non-compact and the compact dimensions of 
our ten-dimensional theory. In principle, this will yield a D = 4 U{N) gauge theory upon 
dimensional reduction. Now, we are again allowed to introduce magnetic fluxes of the form 
(|B.ll ) in the internal dimensions of the D9's worldvolume, again breaking the U{N) gauge 



theory. To each factor Fa we associate a dynamical D9-brane a with gauge group U{l)a- 
Notice that in general a dynamical D9-brane will be composed of ka = rank Fa > 1 type 
IIB D9-branes. In fact, this is how we achieve gauge group rank reduction in this setup. 



^^See the main text for a less abstract definition of slope in terms of magnetized compactifications. 
^^Such gauge unhiggsing occurs if we have in addition A^i = Aa, i = 1, . . . ,Na, i.e., if all the Wilson 
lines are equal. 

''^Recall, however, that our previous discussion is general and describes as well heterotic compactifications 
with non-trivial gauge fields A. 



After turning on F, the spectrum of the dimensionahy reduced theory can be arranged 
in several sectors: 

• Closed strings 

This will provide the gravity sector of the D = 4 theory. 

• Open D9aD9a strings 

Strings beginning and ending on the same dynamical brane will yield a U{l)a gauge 
theory. A stack of dynamical D9-branes with the same flux F^i will yield a gauge 
group U{Na) if there is no relative Wilson lines between them. 

• Open D9aD9^ strings 

Strings stretching between two stacks of N^, dynamical branes may yield D = A 
chiral fermions, as well as scalars, transforming in the bifundamental representation 
{NaiNp) of the corresponding U{Na) x U{Np) gauge group. 

Let us now analyze the supersymmetry conditions above in terms of this physical 
picture. The closed string sector of the theory is mainly sensitive to the geometry of 
the compact manifold Mq. Hence, the condition of being a Calabi-Yau three-fold 
translates into a (tree-level) D = A supersymmetric gravity sector. 

Now let us introduce a magnetic flux satisfying eq.( [B.5| ), that is, let us consider F being 
a (1, l)-from. Without loss of generality, we can consider F to be of the block-diagonal 
form ( [B.ll| ) with a total of K blocks. Each block a corresponds to a stack of Na dynamical 
branes, that is to a flux of the form Ijv^, hence a U{Na) gauge theory. The total flux 
F will provide us with a low energy field theory with gauge group 

Ua=i U{Na). Moreover, 

if each block a satisfies, by itself, the condition ( |B.6| ) we will actually have a low energy 
gauge theory which is a product of U {N^) SYM theories. 

Now, even if that is the case, the D9aD9p spectrum may not be supersymmetric. 
Indeed, each block a will have a slope Cq, = fJ-{Fa). If Ca = cp, then the D9aD9p sector 
will be supersymmetric. To each chiral fermion in a D9aF)9^ sector there will correspond 
a massless complex scalar arising from the same sector, completing a. D = 4 M = 1 chiral 
multiplet in the {Na,N[^) representation. On the other hand, if Cq, ^ Cjs, then the D9aD9i3 
sector will break supersymmetry, and the formerly supersymmetric partners of the chiral 
fermions will be either massive or tachyonic. 

We then see that the conditions for supersymmetry, namely a covariantly constant 
spinor and the two Hermite Yang-Mills equations, can be matched with the supersymmetry 
conditions for the three different sectors of the theory at hand. They can be understood 
from demanding supersymmetry in the gravity, gauge and chiral sectors. 

B.3 Mirror symmetry and string corrections 

Actually, it turns out to be quite instructive to understand these supersymmetry conditions 
in terms of the T-dual or mirror symmetric picture. Indeed, T-duality relates type IIB 
superstring theory on the compact manifold with type IIA theory on its mirror Wq. 
Since we are supposing H = d(l) = in both sides, the only way to achieve a supersymmetric 



gravity sector is to demand that both manifolds admit a covariantly constant spinor, that 
is, that they satisfy the CY3 condition. 




Figure 9: Artist's view of open string mirror symmetry, which relates type IIB magnetized 
compactifications with type IIA intersecting D6-brane worlds. In general, A4e and We will have 
different topology. 

We now consider the open string sector of the theory. The 'mirror' D-branes of a 
magnetized stack of D9-branes will be given by stacks of Na D6-branes wrapping 3-cycles 
Ha on We, and expanding the four non-compact dimensions. The first hermitian Yang- 
Mills condition, eq.( |B.5| ), corresponds in this picture to imposing 11^ to be a Lagrangian 
cycle, that is, to satisfy uj\ua = 0) where u is the Kahler form in Ws. On the other 
hand, the second hermitian Yang-Mills condition ( [B.6D will translate to Ha being a Spe- 
cial Lagrangian (SL) submanifold, which in addition to 11^ being Lagrangian imposes the 
condition lm(e*^"0)|na = 0- Here stands for the holomorphic (3, 0)-form of Wg and 
9a is a constant phase. Finally, the chiral matter in the bifundamental representation 
{Na, Np) arises again in the DGaDQ/s sector of the theory. More precisely it is localized in 
the intersection points of the submanifolds Ha and [37]. See figure ^. 

Again, we find that supersymmetry may be broken or not in different sectors of the 
theory. If we consider Wq being a CY3 and Ha being a Special Lagrangian submanifold 
for every stack a of Na D6-branes, then at low energies we will recover a D = 4 M = 2 
supergravity sector and a gauge sector of U{Na) M = 1 SYM theories. However, if we 
consider the matter content at the intersection of, say, SL's 11^ and 11^ this may yield a 
supersymmetric spectrum or not, depending on the respective phases 9a and O/^. If 9a = 9p 
the DGaDQjS chiral matter will be arranged in AA = 1 supermultiplets, whereas if 9a 7^ 9p 
supersymmetry will be broken in this sector. From the point of view of the effective field 
theory, this can be understood as the appearance of a non- vanishing Fl-term [22,44,45]. 
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In the rest of the paper we will be dealing with one of the most simplest cases of 
open string mirror symmetry. Namely, both and We will be T-dual 6-dimensional tori. 
Moreover, in the type IIB picture we will be considering D9-branes with constant magnetic 
fluxes, which in the type IIA picture correspond to flat D6-branes intersecting at angles. 
This class of configurations have been analyzed in great detail. In particular, in [46] the two 
supersymmetry conditions for D6-brane wrapping a 3-cycle 11 in the intersecting picture 
were understood in terms of F- and D-flatness conditions on the worldvolume gauge theory 
of such D6-brane. This statement should also hold in the D = 4 reduced theory and, by 
mirror symmetry, it should be matched with the hermitian Yang-Mills equations (|B.5| ), 
( B.6 ). We summarize the interpretation of the supersymmetry conditions in Table ||. 



MEW 


IBW 


Field Theory 


Mg = CYs 


We = CY3 


SUGRA 


Fah = 


t^ln, = 


F-flatness 


g'^^F^l = c • 1 


Im(e^^«17)|n, = 


D-flatness 



Table 1: Supersymmetry conditions for both Magnetized Brane Worlds (MBW) and Intersecting 
Brane Worlds (IBW) in terms of the D — A effective field theory. Notice that, in the MBW picture, 
the F-flatness condition imposes a constraint on the complex structure oi Aie, whereas D-flatness 
concerns the Kahler structure. An opposite statements holds for the IBW mirror picture. 

It is important to notice that, although the hermitian Yang-Mills equations and the 
calibration conditions match qualitatively in both sides of the mirror map, there will be a 
quantitative mismatch away from the limit of large volumes and diluted fluxes (or, equiv- 
alently, small angles). Indeed, the notion of /x-stability and the hermitian Yang-Mills 
equations have been derived in the supergravity approximation of string theory. When 
leaving such regime, we would expect stringy a' corrections to the hermitian Yang-Mills 
conditions. Indeed, by considering the conditions for unbroken symmetry from the Dirac- 
Born-Infeld action instead of the SYM Lagrangian, the condition ( [B.7| ) is replaced by the 
MMMS equation [47] 

Ime^^Tr(w + ia'F)" (B.12) 

which involves higher powers of the field strength F. These and other considerations led to a 
modification to the concept of //-stability, namely the Il-stability proposed in [48] . Roughly 
speaking, this proposal amounts to considering the stability condition in the intersecting 
D-brane picture and translate it back to the holomorphic bundle picture by means of the 
mirror symmetry map. As explicitly seen in [14, 16] in toroidal compactifications this 
implies substituting the slope which is related to the tangent of the angles between two 
D-branes, by the angles themselves. Due to this fact, our results concerning three-point 
functions should be only considered as accurate from the string theory point of view only in 
the limit of large compactification volume and diluted fluxes (which corresponds to small 
angles), where the ten-dimensional effective field theory captures all the physics of the 
underlying string theory. 
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